What does inflation really predict? 
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If the inflaton potential has multiple minima, as may be expected in, e.g., the string theory 
"landscape" , inflation predicts a probability distribution for the cosmological parameters describing 
spatial curvature (fitot)j dark energy (pa, w, etc.), the primordial density fluctuations (Q, n s , 
dris/dlnk, etc.) and primordial gravitational waves (r, n t , etc.). We compute this multivariate 
probability distribution for various classes of single-field slow-roll models, exploring its dependence 
on the characteristic inflationary energy scales, the shape of the potential V and and the choice of 
measure underlying the calculation. We find that unless the characteristic scale A(j> on which V 
varies happens to be near the Planck scale, the only aspect of V that matters observationally is 
the statistical distribution of its peaks and troughs. For all energy scales and plausible measures 
considered, we obtain the predictions fitot ~ I ± 10~ J , w = —1 and pA in the observed ballpark 
but uncomfortably high. The high energy limit predicts n 3 « 0.96, dn B /d\nk ~ —0.0006, r ~ 0.15 
and nt ~ —0.02, consistent with observational data and indistinguishable from eternal V oc <f> 2 
inflation. The low-energy limit predicts 5 parameters but prefers larger Q and redder n 3 than 
observed. We discuss the coolness problem, the smoothness problem and the pothole paradox, which 
severely limit the viable class of models and measures. Predictions insensitive to pre-inflationary 
conditions can arise either from eternal inflation attractor behavior or from anthropic selection 
effects probing merely a tiny non-special part of the initial distribution. We argue that these 
two mechanisms are severely challenged by the coolness problem and the smoothness problem, 
respectively. Our findings bode well for detecting an inflationary gravitational wave signature with 
future CMB polarization experiments, with the arguably best-motivated single-field models favoring 
the detectable level r ~ 0.03. 



I. INTRODUCTION 

The spectacular recent progress in observational cos- 
mology has produced measurements of about a dozen cos- 
mological parameters [1-3] as summarized in Table 1 and 
Figure 1 (vertical bands). This observational progress is 
in stark contrast to the theoretical situation: we still lack 
a firm a priori prediction for any one of them. Inflation 
is promising in potentially being able to explain a strik- 
ingly large fraction of these numbers: 8 of the 11, includ- 
ing (as we will argue) p\ and w, the only holdouts be- 
ing the baryon parameter £b (awaiting an understanding 
of baryogenesis) , the dark matter parameter £ c (perhaps 
awaiting an understanding of supersymmetry parameters 
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FIG. 1: Yellow/light grey cosmological parameter distributions 
show inflationary predictions for our §5 example with characteristic 
</>-scale m h ?s 2mpi, F-scale m v = 0.002m P i and uniform measure 
over initial 0-values. Green/dark grey regions show observational 
constraints (ler) from Table 1 and [2, 3]. pa is in Planck units. 



Table 1: Constraints on fundamental cosmological parameters 
computed from WMAP+SDSS as in [2, 3]. Inflation may predict 

the first 8. We have replaced the conventional parameters 
(SIa, Oj,, f^cdmi by the less anthropocentric ones (pa, £b> £c, £k) 
since the latter do not depend on the observing epoch. 
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-1 ± 0.1 




Q 


Scalar fluctuation amplitude S H on horizon 


(2.0 ± 0.2) x 


10~ 5 




Scalar spectral index 


0.98 ± 0.02 






Running of spectral index dn s /d In k 


l«l ^ °- 01 






Tcnsor-to-scalar ratio (Q t /Q) 2 


< 0.36 




n t 


Tensor spectral index 


Unconstraine 


d 


ib 


Baryon mass per photon p h /n^ = jn-p-q 
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governing dark matter freezeout) and the neutrino pa- 
rameter (awaiting an understanding of the origin of 
neutrino masses). 

So what does inflation predict for the parameter vector 

P = (^tot, PA) w,Q,n S) a, r,n 4 )? (1) 

The rhetorical answer is that the answer is determined by 
the inflaton potential V(<p), where is a vector of one or 
more scalar fields, and that the only fairly robust predic- 
tions are f2tot ~ 1 and the slow-roll consistency relation 
r = — 8«t — with known exceptions even to these, e.g. 
[4-7]. Since a bewilderingly large number of potentials 
have been studied in the literature (see [8-15] for recent 
reviews), producing predictions in pretty much ever cor- 
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FIG. 2: A small segment of the inflaton potential V(<j>) for simula- 
tion number 26140213 from §5, with vertical and horizontal energy 
scales mjj = 3m and ra/, = 0.007m, respectively (m = mpi/v / 87r). 
A half-basin of attraction stretches from a local maximum (circle) 
to a local minimum (triangle). The thicker curve indicates regions 
where the slow-roll approximation is valid, with squares indicating 
inflation cndpoints. Many starting points <f> fail to produce galax- 
ies, either because the slow-roll approximation in invalid and there 
is no inflation (three-pointed star), because inflation is rapidly fol- 
lowed by pa < recollapse (five-pointed star) or because inflation 
never ends (six-pointed star). Only a small fraction of starting 
points (like the four-pointed star) give inflation followed by a vac- 
uum density |pa| near zero. 



ner of the currently allowed parameter space, it is easy 
to come away with the impression that inflation is not 
very predictive and hence difficult to test or rule out. 

However, this pessimistic conclusion rests on a rather 
dubious premise: that the way inflation occurred in our 
region of space is the only way that inflation occurred 
anywhere. This is the antithesis of the "landscape" 
paradigm now emerging from string theory, with an ex- 
tremely complicated potential with perhaps I0 300 differ- 
ent minima [16-29], c.f. [30, 31]. 

Whatever the correct fundamental theory may turn 
out to be, one can easily imagine it giving rise to an 
inflaton potential V{<f>) with multiple minima, as in the 
example shown in Figure 2. As long as (p takes on a range 
of values early on , different post-inflationary domains 



will correspond to <p having rolled down into different 
minima, and so the parameter values we measure will 
depend on which domain we inhabit. In other words, 
what is traditionally referred to as the inflaton potential 
is merely the basin of attraction around our particular 
local minimum (fi . Our observed cosmological constant 
will be p\ = V(0 O ), our gravity wave amplitude will be 
given by V(4>\) at the point <p 1 about 55 e-foldings before 
our inflation ended [36, 37], n s will be determined by the 
first and second derivatives of In V at at cf> 1 , and so on. 
Even if the potential has only one minimum, parameters 
such as n s will depend on the direction from which one 
rolls down to this minimum (from the left or from the 
right for single-field inflation like in Figure 2, or from a 
continuum of directions for multi-field inflation where <fi 
is a vector). 

If one or more explicit effective potentials V(<p) emerge 
from string theory or some other fundamental theory [38] , 
then it is clearly worthwhile computing or estimating the 
corresponding 8-dimensional parameter probability dis- 
tribution f p (p) (Figure 1 shows an example), since this 
will provide a powerful observational test of the theory. 
The purpose of the present paper is to lay the groundwork 
for such a calculation in the context of classic slow-roll in- 
flation, computing f p (p) for various classes of randomly 
generated potentials and exploring its dependence both 
on the potential V(0) and on the choice of measure un- 
derlying the calculation. We will see that a number of 
general conclusions can be drawn rather independently 
of the detailed nature of the potential V((j>), depending 
mainly on energy scales. For related discussions of such 
issues in the context in inflationary flow equations, see 
[39-41]. 

We find that the issue of what measure to use when 
sampling 4D spacetime is no less important than the is- 
sue of what potential V(<j>) to use. This measure prob- 
lem cannot be brushed aside as irrelevant philosophy, 
since the measure dramatically affects the testable pre- 
diction for / p (p), and we will argue that some plausible- 
sounding and oft-discussed measures like synchronous 
volume-weighting are already ruled out by observation. 
There is some correct measure that nature subscribes to, 
and we need to figure out which one it is, just as was 
successfully done in the past for the measures allowing 
us to compute probabilities in statistical mechanics and 
quantum physics. 

The rest of this paper is organized as follows. Sec- 
tion II briefly summarizes how the potential V(<p) af- 
fects what we observe, emphasizing that modulo a slight 



1 Having <fi take on a range of values early on is easy to arrange. 
The same quantum fluctuation process responsible for generat- 
ing our observed density fluctuations gives rise to the familiar 



diffusion process that can drive cf> uphill as well as downhill [32] , 
effectively populating all 0-values in a quantum superposition. 
Decoherence [33] then insures that we can for all practical pur- 
poses treat <p as having different classical values in different spa- 
tial regions [34, 35]. Alternatively, an early hot state with en- 
ergy in excess of typical V-values would naturally populate all 
0- values. 
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caveat in the multidimensional case, all inflationally ob- 
servables (including the eight inflation-related parame- 
ters in Table 1) are determined by the cosmic density 
history p(a) alone, independent of the details of V(<p). 
Sections III and IV discuss how the measure affects the 
parameter probability distribution f p (p), i.e., what we 
should expect to observe given V((p), classifying various 
plausible and ruled-out measure candidates. Section V 
presents a suite of numerical experiments exploring the 
effect of inflationary energy scales, and Section VI shows 
how to compute f p (p) from these results when including 
selection effects. Section VII complements our numeri- 
cal exploration with a suite of analytic calculations, and 
Section VIII summarizes our conclusions. Appendix A 
and B give technical details of how we compute p from 
V. 



II. HOW THE POTENTIAL AFFECTS WHAT 
WE OBSERVE 

Given the measure subtleties discussed in the next sec- 
tion, it is crucial to be clear on what question we are 
asking. In this section, we address the following ques- 
tion: If we live in a part of spacetime where the inflaton 
has evolved as <p(a) and the density has evolved as p(a), 
then what do we observe? This is mainly a review of stan- 
dard material, organized to clarify how in many cases, all 
observables can be calculated directly from p(a) alone, 
without needing to know <f>(a) or the inflaton potential. 

We wish to make predictions for four basic observables: 

1 . The curvature of space given by f^tot ■ 

2. The dark energy density as a function of cosmic 
scale factor, p^{a). 

3. The primordial power spectrum of gravitational 
waves 6ff(k). 

4. The primordial power spectrum of adiabatic den- 
sity fluctuations Sn(k). 

It is popular to approximate the three cosmological func- 
tions p\{a), 5jj(k) 2 and Snik) 2 by linear and quadratic 
functions in log-log space, parametrized by parameters 
2-8 from Table 1: 

PA(a) - Pax (a/a Q )-^ 1+w \ (2) 
6 T H (k) 2 = Q 2 xr[k/H) n \ (3) 
8 H (k) 2 = Q 2 x Mfc/ff)/2 ; ( 4 ) 

where H = a/a is the current Hubble parameter, so that 
the wavenumber k = H corresponds to a fluctuation on 
the current horizon scale 23 . Here and throughout, we will 



2 The horizon is c/H — throughout this paper, we use standard 
particle physics units where h = c = k^ = 1, so the Planck mass 
is simply mpi = G~ */ 2 . 

3 For simplicity, we have defined the normalizations (Q, rQ) and 
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FIG. 3: The spatial curvature parameter Qtot, the primordial grav- 
itational wave power spectrum <5|J (k) and the dark energy function 
PA (a) are all determined by the single curve p(a) above, as is the 
primordial density fluctuation spectrum Sff(k) in many cases. Dot- 
ted diagonals of slope —2 are lines of constant comoving horizon size 
a/H -1 = d, governing when fluctuations exit end reenter the hori- 
zon. Inflation (a > 0) is when the cosmic density history p(a) has 
a shallower slope than these diagonals, dlnp/dlna > — 2, and this 
logarithmic slope normally becomes —3 during reheating, —4 dur- 
ing radiation domination, —3 during matter domination and > —2 
now during dark energy domination. When two points lie on the 
same diagonal (like the two triangles), it means that that the hori- 
zon volume at the two epochs is the same comoving spatial region. 
Qtot is constant on these diagonals, approaching unity towards the 
upper right, and unless the p(a) curve crosses the leftmost heavy 
diagonal, we should observe curvature |fitot — 1| > 10 — 5 . The dark 
energy function p\(a) is simply the curve p(a) at late times minus 
the matter contribution. The primordial gravitational wave power 
spectrum is simply the curve p(a) at early times, rescaled vertically 
(by a factor of two since Sjj(k) oc p 1 / 2 ) and horizontally (mapping 
a into k matching horizon exit and reentry with the diagonal lines, 
as with the pairs of triangles, squares and circles). For single-field 
inflation and many multi-field cases, the primordial density fluctu- 
ation spectrum is determined from the curve p(a) and its derivative 
by 8 H (k) oc p 1 / 2 (-dlnp/dlna)- 1 / 2 . 



assume the standard inflationary fluctuation-generation 
mechanism which produces adiabatic fluctuations — for 



tilts (n s ,nt) on the horizon scale. A popular alternative con- 
vention [1—3, 42] is to specify them on the comoving scale cor- 
responding to k = 0.05/Mpc at the present epoch — this corre- 
sponds to shifting the emphasis from, say, scales that are N = 55 
e-foldings before the end of inflation to N ~ 50, making only a 
minor difference in p that is irrelevant for the purposes of this 
paper. 
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a review of alternatives, see [8-15, 43]. 

Figure 3 is an attempt to capture all the inflationary 
observables in one plot. It illustrates how the late-time 
behavior of p(a) encodes the dark energy function Pa (a) 
whereas its early-time behavior determines the curva- 
ture parameter O to t, the gravitational wave power spec- 
trum 5jj{k) and, in many cases, the power spectrum 
of density fluctuation 5n{k). Key epochs in Figure 3 
are a s t a rt (when inflation began; a s t a rt = if there was 
no beginning), a cx it (open triangle, when our current 
Hubble volume left the horizon), (when a fluctua- 
tion of wavenumber k leaves the horizon; open square 
or circle, say), a en d (star, when inflation ended), a re heat 
(cross, when reheating ended), a oq (filled square, matter- 
radiation equality) and ao (filled triangle, now). Figure 3 
also provides a simple graphical way of reading off bounds 
on the number of e-foldings N between the creation of 
our observed horizon-scale fluctuations (open triangle) 
and the end of inflation (star) — most popular models 
give N ~ 55 ±5 [36, 37]. 

For comparison, the ekpyrotic model [44, 45] has fluc- 
tuations leave the horizon during a contracting phase 
when we move up to the left in Figure 3 with steeper 
slope than —2. 

We will now summarize the calculation of the cosmo- 
logical functions pa (a), 5jj{k), 5]j(k) and the eight in- 
flationary parameters from Table 1. As detailed in Ap- 
pendix A, our three cosmological functions and the Hub- 
ble parameter are given in terms of p(a) by 



by (see Appendix A) 



p A (a) = p(a) - /0 7m k(a), 
5 T H {k) = 



8p(ak) 



757r 2 m 4 : 




P(ak) 



757r 2 m 4 (lnp)'' 
din a 



8irp(a) p(a) 



1/2 



dt 



3mpj 



m 



(5) 
(6) 

(7) 

(8) 



where ' = d/dN = —d/d\na, p 7m k(a) is the density con- 
tribution from photons, matter and curvature, the mode 
exit scale is 



afc 



^cxit 



H(a cxit ) p(a cx it) 1/ 



l/2 fc ' 



(9) 



a C xit is the solution to the equation o 2 xit p(a ox i t ) 
agp(ao), and the reduced Planck mass 



m = 



m P1 



mpi 
8tt 5 



(10) 



is defined so as to minimize the number of 87r-factors in 
this paper. The corresponding eight observable cosmo- 
logical parameters in Table 1 are given in terms of p(a) 



|O to t-l| ~ 1(T 5 + 
PA = lim p(a) 



OcxitP(acxit) 
«startp(dstart) 



w 



Q > 



p(a e xit) 



757r 2 m 4 (lnp)'' 



(lnp)' 



r < 8(lnp)', 
n t = -(lnp)', 



(11) 
(12) 
(13) 

(14) 
(15) 

(16) 

(17) 
(18) 



where 



din a 



. The usual slow-roll parameters are 



e = ^(Inp)', 

n \/ ( ln P)" 



in terms of which we have the standard relations 

n s = 1 - 6e + 2?7, 
n t = -2e, 
r < 16e, 



a 



16er7- 24e 2 - 2£ 2 . 



(19) 

(20) 
(21) 



(22) 
(23) 
(24) 
(25) 



As detailed in Appendix A, the inequalities (14), (17) 
and (24) become equalities for single-field inflation and 
also to good approximation for most multi-field mod- 
els. Similarly, the expressions for n s and a are exact 
for single-field models and usually good approximations 
for multi-field models. 

In addition to determining these numerical parameters, 
the p(a)-curve in Figure 3 readily allows qualitative con- 
clusions to be read off: To solve the horizon problem, the 
curve must cross the dotted diagonal that runs through 
the filled triangle, and equation (11) shows that to solve 
the flatness problem, it must also cross the diagonal line 
lying about a factor of 10 further to the left. 



III. HOW THE MEASURE AFFECTS WHAT 
WE OBSERVE 

In this Section, we discuss how to confront an inflation 
model with observational data, i.e., how to compute the 
theoretically predicted probability distribution f p (p) for 
the cosmological parameters once V{<fi) is known. 
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A. Overview of the measure problem 

Inflation helps create an interestingly complex space- 
time, one or more parts of which contain observers 
like us. Loosely speaking, the probability distribution 
fpip) for the cosmological parameters is simply an eight- 
dimensional histogram, reflecting the distribution of pa- 
rameters measured by the different observers. Let us 
now make this more rigorous. To make inflation a fully 
specified and testable theory from which / p (p) can be 
computed, three things need to be specified: 

1. the choice of reference objects (say points, protons, 
planets or observers), 

2. the order in which to count them (this matters 
when there are infinitely many), 

3. the initial conditions (irrelevant in some cases). 

We will refer to the specification of these three things as 
a choice of measure. We will first summarize these three 
issues briefly, then revisit them in greater detail. 




-10 o 

Comoving x- 



1. Reference objects 

The first of the three choices is the most obvious and 
straightforward to deal with, so we will not dwell on it 
here — see [46, 47] for a detailed discussion. The re- 
sulting parameter probability distribution f p (p) depends 
strongly on whether the reference objects are points, pro- 
tons, dark matter halos, galaxies, stars or self-aware ob- 
servers (whatever that means), so the key point is simply 
to be explicit about what reference objects are used and 
to quantify the sensitivity of the results to this choice 
— we will explore various choices below, in Section VI. 
Mathematically, this is simply standard use of condi- 
tional probabilities: we compute the probability distri- 
bution for the parameters p measured from a spacetime 
point given various constraints, say that it is the location 
of a galaxy. 4 



4 Discussions of selection effects often turn heated when some- 
body mentions the "A- word", anthropic. The author feels that 
discussions of the so-called anthropic principle [63-66] have gen- 
erated more heat than light, largely because of a preponderance 
of different and mutually incompatible definitions and interpre- 
tations of what it means. The author is not aware of anybody 
disagreeing with what might be termed MAP, the "minimalistic 
anthropic principle": When testing fundamental theories with 
observational data, ignoring selection effects can give incorrect 
conclusions. This is all we subscribe to in the present paper. 
We wish to test any fundamental theory that predicts V(<f>) by 
calculating its cosmological parameter predictions f p (p), and in- 
cluding selection effects is clearly not optional. The question of 
precisely which selection effects to use (i.e., what reference ob- 
jects to condition on) is difficult and not settled. The appropriate 
response to this is clearly not to give up and ignore selection ef- 
fects altogether, but rather to explore a range of options and 
quantify the extent to which the choice affects the results. 



FIG. 4: Spacetime foliation examples. Solid curves (top and mid- 
dle panels) show infinite hypersurfaces corresponding to the end of 
inflation, starting at t = 0, a = 1 with an inflaton field </>(r) = ex 
with e so small that spatial gradients are negligible. Dashed curve 
(middle panel) corresponds to the present density. The four shaded 
regions of comoving space fail to produce any Hubble volumes 
resembling ours, either because inflation never ends (yellow/light 
grey) or because inflation never starts or lasts less than 55 e-foldings 
(cyan/grey). Lower panel shows the inflaton potential used, with 
vertical dotted lines showing maxima, minima and inflation end- 
points. 

2. Ordering 

Figure 4 illustrates what we will call the ordering prob- 
lem. We have chosen our reference objects (protons, say), 
and need to locate all of them in the spacetime manifold. 
Inflation creates photons and other particles on spatial 
hypersurfaces corresponding to the end of inflation, illus- 
trated by the U-shaped curves in the figure. As detailed 
in Section VII, each of these hypersurfaces typically has 
infinite physical volume, producing infinitely many ref- 
erence objects. Such a hypersurface corresponds to a 
particular basin of attraction in V(<p), i.e., everywhere 
in this comoving volume, the inflaton field <fr rolls down 
to the same local minimum <p . We will usually refer to 
such spatial domains as thermalized regions, pocket uni- 
verses [48-50] or "pockets" for brevity — other terms for 
them include bubble universes and Level I multiverses 
[51]. 5 If we order the reference objects 1, 2, 3, ... and let 



J Note that there is not a one-to-one correspondence between 
basins and pockets. Each pocket corresponds to exactly one 
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Pi denote the parameters observed from the i object, 
then we can compute the probability distribution as 



/ p (p)= lim - V^(p-Pi). (26) 

n — >oo n. * — ' 



The crux is that when there are infinitely many objects 
[48-50], the ordering affects the answer! For example, 
suppose we order the objects by increasing physical dis- 
tance from some given point r at fixed time (defined as 
fixed cosmic density). This may sound like an innocent 
and reasonable algorithm, since it is equivalent to the fa- 
miliar procedure of averaging in a spherical volume and 
then letting the sphere radius approach infinity. How- 
ever, since the pockets (thermalized regions) arc infinite, 
this ordering will only include objects within the same 
pocket that contains the reference point r, giving zero 
statistical weight to the reference objects in all other 
pockets. 

This illustrates a difficult aspect of the ordering prob- 
lem: how to decide the relative order of reference objects 
in different pockets. This is equivalent to assigning sta- 
tistical weights to the different pockets once the ordering 
and hence the p-distribution has been determined within 
each pocket. As detailed below, many attempts to solve 
the ordering problem have implicitly involved grouping 
the reference objects by the value of some physical pa- 
rameter at their spacetime location [52-61], say by proper 
time t (Figure 4, top panel), scale factor a (Figure 4, mid- 
dle panel) or density p (roughly the black curves in top 
two panels). 

In Section IV, we will review and discuss various or- 
dering choices and their implications. This is a severe 
problem in the sense that the ordering strongly affects 
the testable prediction / p (p), yet we currently lack a 
compelling argument leading to a unique choice. As de- 
tailed below, it even affects the sense in which we should 
think of inflation as eternal. This embarrassing problem 
is clearly one of physics rather than philosophy, since its 
resolution affects the testable inflationary predictions for 
/ p (p). On an optimistic note, the measure problem (how 
to compute probabilities) plagued both both statistical 
mechanics and quantum physics early on, so there is real 
hope that inflation too can overcome its birth pains and 
become a testable theory whose probability predictions 
are unique. 



basin, but a basin can correspond to multiple pockets, either 
because of initial conditions (in Figure 4, imagine, say, an initial 
inflaton field <^(r) = <f> pca k + esina; that is draped over the max- 
imum at <^p Ca k multiple times) or because of quantum diffusion 
(which can drape the inflaton field repeatedly over the maximum 
even starting with the simple initial condition </>(r) = ex where 
the field crosses the maximum only once) [32]). 



3. Initial conditions and predictions 

Let us cast the above issues in mathematical form, in- 
cluding the initial conditions <fr which are the inflaton 
field and its first derivative early on: 



Given any particular initial conditions 3? and inflaton 
potential V(<p), the standard inflation formalism allows 
us to compute the probability distribution / p (p; <J>, V) 
for the parameters p. If quantum diffusion is negligible, 
then the parameters are uniquely determined by $ and 
V and we have simply 

/ p (p;*,V0w*tP-P(*.^)]. (28) 

where p(3>, V) is computed as in Section II or by go- 
ing beyond the slow-roll approximation if necessary [11]. 
When quantum diffusion is important, a given starting 
point $ no longer gives a unique classical time-evolution 
4>(t) and so the probability distribution f p (p;$,V) 
widens — for this case, / p (p; «&, V) can be evaluated ei- 
ther by solving a stochastic ordinary differential equation 
or, equivalently, by integrating a Fokkcr-Planck equation 
[8, 53, 62, 103-106]. As will be discussed in Section IV, 
there are even some choices of measure where attractor 
dynamics [8, 53, 62, 103-106] causes / p (p;*,V) to be 
completely independent of the initial conditions 

Quite generally, we can now write the theoretical pre- 
diction for the parameter probability distribution as 

/p(p)=Kp)«/ P (p;*,V0>3>v, ( 29 ) 

where w(p) is a weight function reflecting the choice and 
ordering of reference objects. Here ( )<j> denotes aver- 
aging over initial conditions <fr, to cover cases where the 
theory predicts a distribution of initial conditions rather 
than some particular initial conditions. Analogously, ( )y 
denotes averaging over inflaton potentials V to cover 
cases where this is necessary — on could, for instance, 
envision theories predicting a quantum superposition of 
many different effective inflaton potentials corresponding 
to different false vacua of the underlying theory. 

For notational convenience, we augment our parameter 
vector from equation (1) to 

p = (Cl tot ,pA,w,Q,n s ,a,r,n t ,N tot ), (30) 

making it include also N tot , the total number of e- foldings 
during inflation. Although, unlike the other parameters, 
iVtot is not directly observable, it enters in many inter- 
esting choices of the weight function w(p). In Section VI 
we will explore a variety of weight functions w(p) includ- 
ing factors such as e 3Ntot (weighting by the thermalized 
volume produced), 0(N tot — 55) (a Heaviside step func- 
tion giving weight only to regions inflating by at least 
55 e-foldings, a requirement linked to galaxy formation), 



7 



/ halo (Q, p A , fitot) (the fraction of protons ending up in 
dark matter halos) and /q(Q) (a factor related to galaxy 
formation and planetary stability). 

We will find that the choice of measure makes a dra- 
matic difference, and leads to one of three qualitatively 
different situations for the theoretical prediction / p (p): 

1. Predictions independent of initial conditions be- 
cause of attractor dynamics. 

2. Predictions independent of initial conditions be- 
cause selection effects in w(p) probe only tiny non- 
special part of initial distribution. 

3. Predictions dependent on initial conditions, i.e., on 
pre-inflationary physics. 

Case 1 is discussed in Section VII, with the conclusion 
that some implementations of it are observationally ruled 
out. Case 2 is discussed in Section V and is a multidimen- 
sional analogy of the classic example of anthropic con- 
straints on pa [61, 64, 67-71] 6 . For some cases explored in 
Section V, we will find the predictions to be independent 
of the details of the inflaton potential V(cf>) as well. Cases 
1 and 2 are nice in the sense that we can make quanti- 
tative inflationary predictions from V(</>) alone, without 
having a theory of pre-inflationary physics. On the other 
hand, case 3 offers an opportunity to learn about pre- 
inflationary physics from cosmological observations of p. 

In summary, we have seen that the when computing 
inflationary predictions, the choice of measure discussed 
in this section has just as important observational con- 
sequences as the choice of inflaton potential discussed in 
Section II. Moreover, we found that the measure prob- 
lem splits into three parts: one straightforward (chosing 
reference objects), one hard (chosing an ordering) and 
one perhaps irrelevant (chosing initial conditions). 



IV. POSSIBLE SOLUTIONS TO THE 
ORDERING PROBLEM 

Above we saw that to make inflation a testable physi- 
cal theory, the ordering problem must be solved. This is 
a wide open problem in dire need of further work — the 
purpose of this section is not to give a solution, merely 
to describe some possible approaches to solving it. We 
will compare their strengths and shortcomings in the re- 
mainder of the paper, first using numerical computations 



Weinberg's argument [68] went as follows: Since galaxies only 
form for \p\\ < 10 , all models predicting an a priori proba- 
bility distribution /(pa) arc equivalent as long as the character- 
istic width of the distribution is 3> 10~ 123 and the value Pa = 
is not in any way special. The probability distribution for p\ 
seen from a reference object (in this case a galaxy) can therefore 
be computed independently of the detailed shape of the initial 
condition distribution f(p\), simply treating /(pa) as constant. 



in Section V and then using analytic calculations in Sec- 
tion VII. We will see that all orderings can be viewed as 
regularization techniques to deal with infinities, and that 
many are observationally ruled out. 

A variety of solutions to the ordering problem have 
been discussed in the literature, albeit using different 
terminology. We will consider two broad classes of or- 
derings: 

1. Global orderings based on some time variable (say 
t or a). 

2. Pocket-based orderings (ordering separately within 
each thermalized region (pocket), then averaging 
the pocket results with some weighting (say equally, 
by 0- volume or by r- volume). 

A. Global time-based orderings 

The first class of orderings foliates spacetime into a se- 
quence of three-dimensional spatial hypersurfaces, each 
corresponding to a fixed "time", and computes the pa- 
rameter probability distribution f P (p) separately on each 
hypersurface. We will see that this is in many cases 
equivalent to simply ordering the reference objects by in- 
creasing formation "time" . The utility of this approach 
stems from two facts: 

• In many cases, the ordering problem vanishes on 
each individual hypersurface. With time variables 
like t or a, a finite initial volume will give a finite 
hypersurface volume for any fixed future time, so 
that all orderings of the (finitely many) reference 
objects give the exact same answer f p (p). 

• In many cases, / p (p) converges to a well-defined 
distribution as this "time" approaches infinity, 
eliminating the need to select a particular hyper- 
surface. 

As detailed below, this first class of orderings is inti- 
mately linked to the infamous problem of chosing a time 
variable in the inflationary Fokkcr-Planck equation, and 
suffers from serious problems [48-50, 52-61]. 

For all cases where the total number of reference ob- 
jects grows exponentially over time, any asymptotic t — > 
oo distribution for reference objects at fixed time will be 
identical to the asymptotic distribution for reference ob- 
jects formed at time < t, because older objects which 
formed at time <C t are exponentially underrepresented. 
The latter distribution is the one corresponding to simply 
ordering all reference objects by their formation time, so 
we will use the expressions "^-ordering" and "f-foliation" 
interchangeably. 

B. Pocket- based orderings 

In contrast, attempting to use density, temperature or 
expansion rate (or indeed any quantity which, unlike t or 
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a, is physically observable) as a time variable to foliate by 
leads to the second class of orderings — such "time vari- 
ables" do not specify the ordering sufficiently for inflation 
to predict a unique parameter probability distribution, 
since there will typically be infinitely many reference ob- 
jects at any fixed "time" . Specifically, infinite numbers 
of reference objects at a given density, say, can typically 
be found in each of multiple thermalized pockets (the 
U-shaped regions in Figure 4, say) and neither the intra- 
pocket ordering nor the inter-pocket weighting has been 
specified — we will discuss multiple options for both. 

1. Intra-pocket orderings 

For the former, i.e., within a given pocket, all orderings 
clearly give the same answer if the volume is finite. We 
will study the following two orderings for the infinite- 
volume case: 

1. Empirical ordering: Order the reference objects by 
increasing physical distance from some given point 
r. 

2. Volume-based ordering: Order the reference objects 
sampling each spatial region in proportion to its 
relative volume fraction, defined in terms of how 
much it has inflated. 

The empirical ordering is equivalent to the procedure 
that an empirically minded observer at r would adopt: 
compute the parameter probability distribution by aver- 
aging in a spherical volume centered at r (on a hyper- 
surface defined by a fixed density, say) and then let the 
sphere radius approach infinity. This procedure has been 
advocated by [72, 73], who applied it on the hypersur- 
face corresponding to the end of inflation. It is easy to 
see that its predictions depend only on which pocket the 
point r lies in, not on its location within that pocket. 

In contrast, a volume-based ordering would sample 
each spatial region with a frequency proportional to how 
much is has expanded since some early reference instant 
— the latter is challenging to define unambiguously, and 
currently there is no rigorously defined implementation of 
a volume-based ordering. These two orderings are distin- 
guished by their answer to the question "are some vol- 
umes more infinite than others" , to which they answer 
"no" and "yes", respectively. Consider, for example, the 
spatial region with < x < 22 in Figure 4. This is a 
single basin of attraction with infinite thermalized vol- 
ume both to the left and right side of the plane x = 14, 
corresponding to the inflaton having rolled down to the 
minimum from the left and from the right, respectively. 
This predicts a cosmological parameter distribution of 
the form 

f P (p) = PlS(p - Pi ) + (1 - PL )S(p - pa), (31) 

where pl and pa are the cosmological parameter vec- 
tors corresponding to (j> rolling down from the left and 



right sides, respectively, and ph is the fraction of the 
reference objects in regions where 4> rolled down from 
the left. The vectors pl and p^ are determined by the 
potential V(<p) alone as per Section II, whereas the con- 
stant pl depends also on the measure. As long as both 
half-basins have infinite volume, the empirical ordering 
will predict pl = 0.5, whereas a volume-based ordering 
will generally not. For instance, if the reference instant 
corresponds to a given density, then volume-based or- 
dering assigns a larger probability to the side where it 
takes more e-foldings to roll down. For the multi-field 
case where is a d-dimcnsional vector, p^ and p# are 
generically replaced by a continuum of parameter vectors 
corresponding to the directions in </>-space from which the 
inflaton rolled down, and equation (31) gets generalized 
to a weighted angular average in 0-space. 

In summary, both the empirical and volume-based 
intra-basin measures agree within any finite sub-volume 
(where they both weight by volume, i.e., reduce to stan- 
dard Lebesgue measure), but differ in how to resolve the 
ambiguities associated with infinite volumes. 

2. Inter-pocket weightings 

Once an intra-pocket ordering has been prescribed 
within each pocket, the inter-pocket weighting must spec- 
ify how reference objects from the different pockets are 
to be merged into a single ordered sequence. For in- 
stance, ordering objects A\, A 2 , ... from pocket A and ob- 
jects Bi, B 2 , ... from pocket B jointly as A\, B\, A 2 , B 2 , ... 
corresponds to giving equal weight to the two pockets, 
whereas A\, A 2 , B\, A3, A 4 , B 2 , A 5 , A 6 , B 3 ... gives twice 
as much weight to pocket A as to pocket B. There are at 
least five obvious candidates for the inter-pocket weight- 
ing: 

f . Equal weighting for each basin of attraction (in <fi- 
space) 

2. Equal weighting for each pocket (in physical space) 

3. Weighting by relative <p- volume 

4. Weighting by relative thermalized physical volume 

5. Weighting by relative number of reference objects 

1 and 2 generally differ, since many disconnected pockets 
can have <p roll down into the same basin of attraction. 
An explicit definition of 4 is given in [74] for the special 
case of inflaton potentials with a single peak. 

Since inflation typically produces infinitely many pock- 
ets, there is a potential subtlety related to pocket order- 
ing within the weighted average. We define two pock- 
ets to be of the same type if they are observationally 
indistinguishable, i.e., give the same / p (p). If there 
are only finitely many types of pockets, the inter-pocket 
weighting is fully specified by giving the weight for each 
type. If there are infinitely many types, however, a naive 
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weighted average over the types is ambiguous, since the 
order in which one adds the types (or pockets) matters. 
This is why, most generally, one needs to merge the ref- 
erence objects from all pockets into a single ordered se- 
quence as stipulated above. 



C. Symmetry arguments 

In the following sections, we will explore the predic- 
tions of various orderings listed above, and find that they 
are sufficiently different that it might be possible to deter- 
mine observationally which one, if any, is the correct one. 
In the spirit of attempting to predict rather than postdict 
the correct answer, we note that there are three heuristic 
arguments favoring orderings where different pockets are 
treated symmetrically, i.e., receive equal weight. 

The first argument is that, mathematically, two sets 
contain the same number of objects if there exists a one- 
to-one correspondence between them. Applying this cri- 
terion to the case of a one-dimensional inflaton poten- 
tial, all countable infinities are equal and hence all in- 
finite half-pockets should get equal weight, which cor- 
responds to empirical intra-pocket weighting and equal 
inter-pocket weighting. 

The second argument involves the "pothole paradox" 
that will be presented in Section VII F, and also favors 
equal intra-pocket weighting. 

The third argument involves the following Gedanken 
experiment. Consider two basins of attraction in a one- 
dimensional inflaton potential V ((/)), and let A, B, C and 
D denote four infinite spatial regions where (p has rolled 
down from the left into basin 1, from the right into basin 

1, from the left into basin 2 and from the right into basin 

2, respectively. In the special case where V(<p) has iden- 
tical shape in the left parts of the two basins, i.e., where 
there is some offset A<p such that V(4>) = V(<p + A</>) 
for all 4> in the first half-basin, then one can argue that 
volumes A and C should receive equal statistical weight 
by symmetry. 7 Since empirical intra-basin weighting 
equates the weight of A with B and the weight of C 
with D, transitivity of equivalence implies that B and 
D have equal weight, and therefore that the two pockets 
A+B and C+D have equal weight even though they are 
causally disconnected. Since this argument applies re- 
gardless of the potential shape in half-basins C and D, it 
suggests that all half-pockets and hence all pockets have 
equal weight. The assumed shape-equivalence of A and 
C can be eliminated from the argument if we consider a 
messy random potential like in Section V with infinitely 
many basins where for any given half-basin, there will 
be another one that approximates its shape arbitrarily 



This argument fails for measures where quantum diffu- 
sion/tunneling from nearby basins is important. 



closely. In summary, this suggests that empirical intra- 
basin weighting implies equal inter-pocket weighting. 

We close this section by noting that the ordering prob- 
lem is closely linked (but not equivalent) to the condi- 
tional probability calculation corresponding to condition- 
ing on a class of reference object. For example, chosing 
our reference objects to be "protons at time i" implies us- 
ing the above-mentioned f-foliation. In contrast, chosing 
reference objects to be "protons at density p" (or con- 
ditioned on T, H or another observable) does not fully 
specify an ordering for the reasons given above, merely 
implying a pocket-based ordering. 

V. MONTE CARLO ANALYSIS 

Above we discussed how the cosmological parameter 
probability distribution f p (p) predicted by inflation de- 
pends on the inflaton potential V(cp) (Section II) and on 
the measure (Section III). We wish to keep our discus- 
sion from becoming overly abstract, so let us now com- 
pute the predictions for some concrete examples. We will 
then quantify the measure dependence for a wider class 
of orderings below in Section VII. 

A. Choice of potentials 

Since the inflaton cp and the potential V have units of 
mass and density (mass 4 ), respectively, let us write the 
potential as 

V(cp) = m v 4 f (*-) , (32) 

where / is a dimensionless function with values and 
derivatives of order unity, and the characteristic energy 
scales of the horizontal and vertical axes in a plot of V(cp) 
are absorbed into the two constants rrih and m v . If V(cp) 
emerges from some fundamental physical theory without 
fine-tuning, one might naturally expect ~ m v , but 
we will also explore the more general case where the hor- 
izontal and vertical energy scales are different. 

Ideally, we would wish to explore generic potentials 
whose shape / is derived from some fundamental the- 
ory. Currently lacking this, we will instead explore po- 
tentials that are generic in the sense of being randomly 
generated. We will explore the case where / is a one- 
dimensional Gaussian random field with unit variance 
and power spectrum P(q) oc q J e~ q I 2 (see Figure 2). In 
practice, we define 




where the (2n + 1) Fourier coefficients a_„, a n are in- 
dependent real-valued Gaussian random variables with 
zero mean and variance 



(a2>=AXy e -V (34) 
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q = k/\fn and the normalization constant A is chosen 

so that (f(x) 2 ) = ±Efc=-n( a fc) = L We take n = 100 
and 7 = as our baseline, but explore a variety of other 
choices below. We will find that the resulting parameter 
probability distribution / p (p) is rather insensitive to the 
detailed shape of /, but depends strongly on the energy 
scales mh and m v . A key goal of this section is to identify 
which features of a complicated potential V(4>) are most 
important in determining the testable predictions / P (p). 

B. Choice of measures 

We will compute results for two different measures that 
we term Measure A and Measure B. Both weight all half- 
basins roughly equally, the difference being that A counts 
only those producing an infinite volume. Measure A thus 
corresponds roughly to the empirical intra-pocket order- 
ing and the equal weighting of basins. As described be- 
low, this also corresponds approximately to a uniform 
distribution over the initial conditions (4>,<p), including 
the volume weighting factor e 3Ntot into the w(p)-factor 
of equation (29) for Measure A but excluding it for Mea- 
sure B. Both of these two measures may be the wrong 
choice, but they are least free from the most egregious 
paradoxes and observational contradictions that afflict 
many other measures as discussed below in Section VII, 
and hence instructive to explore quantitatively. 

Specifically, we perform the following numerical exper- 
iment: 

1. Generate a random potential V(4>), start at <f> = 0. 

2. Go uphill until a maximum is reached. 

3. If V < (Big Crunch imminent) or the slow-roll 
approximation (SRA) is invalid (e > 1 or |7/| > 1), 
count as an Ar tot = failure and stop, otherwise 
evolve according to the SRA. 

4. If (j) gets stuck in an eternally inflating local mini- 
mum, count as a failure and stop (will receive zero 
weight since no reference objects produced — not 
even protons). 

5. When the SRA breaks down (which defines a en d), 
compute p using equations (11)-(18), taking a cx i t 
to be 55 e-foldings earlier (a ex it = e _55 a on( j) and 
evaluating p\ — V(<p) at the subsequent local min- 
imum of the potential. 

6. Repeat to obtain a statistically large sample, say 
10 9 times. 

7. Repeat the entire experiment for a range of energy 
scales mh and power spectrum shape parameters 7. 

Step 2 is optional: we include it for Measure A but not for 
Measure B. We are thus simulating "shotgun inflation" in 
the sense that we spray starting points randomly across 
the potential surface (Measure B) or at random maxima 



(Measure A). In Section VI, we will use either galaxies 
or objects in galaxies as reference objects, so the weight 
function w(p) from equation (29) will eliminate all cases 
where there is too much dark energy (\p\\ 3> 10~ 112 Q 3 ) 
for galaxies to form or where black holes prevent the for- 
mation of reference objects (for V to t 55 as per Sec- 
tion VIC 2). As detailed in Section VII, Measure A will 
automatically give N tot = 00, whereas Measure B will 
give Nt t — 00 only if the maximum above the starting 
point satisfies the SRA, i.e., has r\ > —1. This V tot - 
cutoff thus affects only Measure B. 

The two measures defined by this prescription were 
chosen to maximize computational efficiency, not to be 
maximally simple to interpret. However, we will now see 
that, modulo some weight factors or order unity, they 
correspond to weighting all half-basins equally and also 
to uniform initial conditions. 

Since the statistical properties of our Gaussian ran- 
dom field V((f>) are translationally invariant, our choice 
to start at <j) = is equivalent to starting at a random 
</>-value drawn from a uniform distribution. Apart from 
the numerical practicality that n is finite so that V is 
periodic, we can thus reinterpret our generation of many 
random potentials as simply trying many different start- 
ing points in the same infinitely extended potential. The 
probability that we roll down into any particular half- 
basin is therefore proportional to its length for Measure 
A and to the length of the sub- region giving 7V t ot > 55 for 
Measure B. The length of a half-basin (the distance be- 
tween a maximum and a minimum of V) is of order mt 
with random variations of order unity, so the different 
half-basins get weighted roughly equally. 

The parameter p\ of course depends only on which 
basin we roll down into, since it is the height of its min- 
imum. Since the potential is one-dimensional, the pa- 
rameters (Q,n s ,a,r,n t ) depend only on the half-basin 
(whether we roll down from the left or from the right), 
so the only parameters affected by where in the basin we 
start are O t ot and N to t — for Measure A, where we start 
at a maximum, they are f2tot = 1 and N to t — 00 • 

Finally, let us discuss the effect of requiring the slow- 
roll approximation (SRA) to hold. When instead using 
the exact inflationary dynamics of equation (A14), we 
find that almost all cases giving substantial amounts of 
inflation do so while the slow-roll approximation (SRA) 
is valid. Because our w(p) gives negligible weight to cases 
with < 55 e-foldings (black holes prevent the formation of 
reference objects as per Section VI C 2), our statistics will 
be dominated by trajectories <p(t) involving at least one 
slow-roll period. Our prescription counts those starting 
with slow-roll and discards those where the SRA becomes 
valid only later on (for instance, a large 4> can be finely 
tuned to just barely make </> roll up a hill, so that the 
SRA becomes valid near the top). If e < 1 and |7j| < 1 at 
a given position </>, then the well-known SRA attractor 
behavior [11] will soon drive the derivative <p to the SRA 
value <f> i=a —ffiV '(</>) y/3/V(4>) if it starts out in the slow- 
roll range \<j)\ < y/2V{4>). If we began with random and 
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!g a / a end 

FIG. 5: The cosmic density history p(a) is shown in the three pan- 
els with labeling on the right side for a handful of simulations with 
Measure A and inflation mass scales (m/,,m„)=(0.2m, 10 _6 m), 
(2m, 0.0004m) and (100m, 0.02m). The corresponding scalar and 
tensor power spectra are shown with labeling on the left side. The 
shaded vertical band indicates the range 43-55 e-foldings from the 
end of inflation where there is current hope to measure these pri- 
mordial fluctuations [75, 76]. 

uniformly distributed initial conditions $ = (</>,</>), the 
fraction of the models surviving the (f> cut would thus 
be ex V 1 / 2 , so that the 0-cut is equivalent to starting 
with the SRA value of </> and weighting the initial (f>- 
distribution by V{<j)) 1/2 . We find that V 1/2 does not 
vary strongly between successfully inflating models, so 
our prescription is roughly equivalent to simply starting 
with uniformly random initial conditions 4>. 




!g a / a end 

FIG. 6: Same as Figure 5, but for Measure B, which includes 
thermalized regions with finite volume. 



For computational efficiency, we ignore quantum dif- 
fusion in our calculations. We will discuss the effects of 
diffusion in Section V G and find that they are unim- 
portant for our qualitative conclusions. As discussed in 
Section VII, diffusion can, with some measures, drasti- 
cally alter the relative weighting between different basins. 
Once we have specified the inter-pocket weighting, how- 
ever, the only observable effect that diffusion can have is 
to modify the evolution during the last 55 or so e-foldings 
— this is important only for models with Q > 1 and is 
completely irrelevant for models with Q ~ 10~ 5 . 
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FIG. 7: The cosmological parameter distributions predicted by 
three inflation models with m v ~ 0.007m,; and Measure A. 
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FIG. 8: The cosmological parameter distributions predicted by 
three inflation models with m v ~ 0.007rrii, and Measure B. 



C. Basic results 



Figures 5 and 6 shows the resulting p(a)-curves for a 
small sample of models for the two measures. From these 
curves, the cosmological functions <$jj(fc) and Sjj(k) are 
readily obtained using equations (6) and (7) — these two 
figures show examples. Our cosmological parameters p 
are computed from the p(a)-curves using equations (11)- 
(18), giving histograms such as those shown in Figures 7 
and 8. Let us now explore how these results depend on 
the inflation mass scales and ra„ and on the param- 
eter 7 controlling the power spectrum of the inflaton po- 
tential. 

The mu-dependence is readily computed analytically. 
Since quantum diffusion is negligible for the cases we 
consider, the vertical inflation scale m v affects none of 
the parameters except the overall amplitudes Q and 
Qt — tQ, which both scale as m v 2 (the pA-distribution is 
affected only indirectly via Q, as will become clear in Sec- 



tion VI). The effect of changing the horizontal scale rat 
(the scale on which the inflaton potential wiggles), how- 
ever, needs to be computed numerically, and is complex 
and interesting for most parameters. We compute the 
parameter distribution f p (p) for the values m^/m =0.2, 
0.5, 1, 2, 3, 5, 7, 10, 20, 30, 100 and 1000, and the results 
are summarized in Figures 9 and 10. For our baseline 
calculations, we take 7 = 0. We repeat the entire grid of 
calculations for 7 = 2 and 7 = 4 and find that this has 
little effect on the results — the analytic results below 
will clarify why. 

Wc find that there are three regimes of horizontal mass 
scale giving qualitatively different behavior: low energy 
to/i <C m, intermediate energy nih ~ fh and high energy 
mh 3> fh. We will now discuss these three regimes in 
turn. 
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FIG. 9: Each shaded band shows the central 68% of a parameter 
probability distribution (such as those in Figure 7) as a function 
of horizontal inflation mass scale m^, for Measure A. All labeled 
bands have the vertical scale fixed at m v = 0.007m. The second 
pair of bands for Q and Qt (unlabeled, above labeled bands with 
same colors/shades) have m v varying so that m v = rrih- 
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FIG. 10: Same as Figure 9, but for Measure B. which includes 
thermalized regions with finite volume. 



D. The high energy limit m/, S> m 

Although we lack a theory of quantum gravity, it is 
commonly assumed that quantum gravity corrections will 
make contributions to the inflaton potential in the form 
of a power series in (c/>/mpi) with coefficients of order 
unity, thus causing V(<p) to vary substantially when on 
a scale Acj) ~ mp\. In other words, one assumes that 
quantum gravity predicts rrih ^ mpi ~ 5m, rendering the 
high-energy regime rrih 3> fh physically unnatural and 
poorly motivated. This classic argument [8-15] is also 
the standard objection to large-field inflation models like 
V((j)) oc (j> 2 , which require this sort of extreme flatness to 
agree with observation — see [77] for a counterarguments 
and discussion of these issues. 

These issues aside, what are the testable predictions in 
the rrih S> fh limit? Figures 9 and 10 reveal an interest- 
ing limiting behavior: when rrih increases far above the 
Planck scale, the predictions for n s , r and a become very 
sharp, with almost no random scatter, independently of 
whether Measure A or Measure B is used. This behavior 
can be traced back to Figures 5 and 6, illustrating that 
as rrih ~ > oo, all the p(a)-cmves have the same shape, dif- 
fering merely in amplitude. The slow-roll parameters e, 
•q and £2, which depend only on the shape of p(a), there- 
fore approach constants and the corresponding probabil- 
ity distributions for r, Tit, n s and a in Figures 7 and 8 
approach Dirac ^-functions as rrih — * 00. 

The reason for this limiting behavior is easy to under- 
stand. Increasing rrih makes it easier to satisfy the slow- 
roll conditions e -C 1 and \r]\ <C 1, since equation (32) 
implies that 
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(36) 



at generic points </>, so that as rrih ~ ¥ 00 > inflation oc- 
curs at almost all </>-values where V((f>) > 0. If inflation 
ever ends, it does so because / drops very close to zero, 
either because of e when / = (m; l /m)|/'|/2 or because 
of rj when / = (mh/fh) 2 \f"\. Recall that the dimension- 
less inflaton potential / is a function whose values and 
derivatives are generically of order unity. In either of 
the two cases, our rrih 3> fh limit therefore implies that 
< / <C 1 when inflation ends, so that the observational 
consequences depend only on the behavior of / over a 
tiny range very close to zero. Since this range is much 
smaller than the natural scale on which / varies, the only 
aspects of / that have any observational consequences are 
the first terms of its Taylor expansion in this range. 

The most likely behavior in the region where / ~ 
is that / is roughly a straight line with slope of order 
unity. For most of our Monte Carlo simulations, infla- 
tion is therefore ended by e = 1 when / ~ fh/rrih, after 
which the inflaton fast-rolls [6] down to negative values 
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of the potential and space promptly recollapses without 
producing any galaxies or other reference objects. 

Although these are the most common cases, they are 
given zero statistical weight by the w(p)-factor in equa- 
tion (29) and hence do not contribute to the predicted 
parameter probability distribution / p (p). The only cases 
that produce galaxies are those where /' is unusually 
small in this region, so that / flattens out and takes 
a local minimum where p\ = V((j>) = m v A f {<f> / mh) is 
close enough to zero to produce galaxies. Since this con- 
straint is roughly |pa| 3> 10~ 112 <2 3 (see Section VI), the 
only segments of a potential V(<f>) that contribute to 
the parameter probability distribution are those whose 
minimum V(4> s t op ) is much closer to zero than the V- 
value where inflation ends, in other words those segments 
that are well-approximated by a parabola V(<f>) oc <ft 2 . 
For this well-studied case [11, 32], we have e = r\ = 
2m 2 /((f> - </> st0 p) 2 , which gives cnd = </> stop ± \f2fh and 
the following cosmological functions and parameters (see 
Appendix B and, e.g., [11]): 



p(a) 

Q 
Qt 



p(a end ){l + 2N), 
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75tt 
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(37) 
(38) 



V2N + 1Q» w 1.5Q* 



a = 



75tt 

4 2 
— — - w l - _ w 0.963, 
2N+1 N 
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0.0006, (41) 



(39) 
(40) 
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n t = - 



2N + 1 N 
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2NTl ~ ~N 



N 2 
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0.02, 



(42) 
(43) 



where N = a end /a as usual, 0* = — \/ fZ;„ and 
/min ~ 1 is the second derivative of the dimensionless 
inflaton potential / at the local minimum into which it 
rolls down. The number of c-foldings N corresponding 
to our current horizon volume is given by [36] 



jV«63.3-— ln^, 

12 Prch 



(44) 



where the last term is 0, 4, 11 and 15 for reheating in- 
stantaneously, at < 10 11 GeV (to avoid overproduction 
of gravitinos in SUSY models), at the electroweak scale 
100 GeV and at the BBN scale 1 MeV (at lower reheat 
energies, no protons are produced — but see [78]). Wc 
have used N = 55 in all our plots. Finally, we expect the 
observed O t ot ~ 10~ 5 as rrih — ► oo, since the decreased 
slow-roll speed increases the total number of e-foldings 
needed to roll down from a given starting point (j>, giving 

Abeforc = ln(a sta rt/aoxit) * OO as fflfr ► oo. 

The above analytic arguments are confirmed by our 
numerical calculations. Figures 5 and 6 (bottom) show 



that the density histories p{a) have the shape predicted 
by equation (37) — the scatter in amplitudes is due to the 
fact that pend oc Q 2 oc / min , and this second derivative of 
order unity differs from one minimum to another. Fig- 
ures 9 and 10 show that the dispersions in n s , r, n t , a, e, 
n and £ 2 drop towards zero as rrih increases, and that 
the constant values they approach are precisely those 
given by equations (40)-(43), i.e., (n s , r, nt, a, e, r], £2) ~ 
(0.963, 0.15, -0.02, -0.0006, 0.009, 0.009, 0). 

Also, Figure 10 shows that the success rate drops 
oc m,h~ 2 as rrih — > 00, as predicted by our above dis- 
cussion: inflation ends with 77 ~ 1 when / = / enc j ~ 
(mh/fh) 2 f" ~ (m/j/m) 2 -c 1 according to equation (36), 
and the probability that the next minimum lies between 
this tiny positive value and zero is thus proportional to 
/end <£m h - 2 . 

Equation (38) shows that reproducing the observed 
fluctuation level Q w 2 x 10~ 5 requires rrih ~ 
0.003^rhm v , so in addition to the possible quantum grav- 
ity problems with having rrih 1, these models also re- 
quire tuning in the sense that m v -C m h (they are far 
from the band marked "natural" in Figure 11). 

In summary, the parabolic inflaton potential V((/>) oc 
4> 2 has been widely studied in the literature, and is a clas- 
sic example giving eternal inflation with an unbounded 
potential [32] . It is also the only integer power law poten- 
tial that is still consistent with observational data now 
that the c/> 4 model is firmly ruled out [3]. In this subsec- 
tion, we have seen that the predictions of this particular 
model hold also for a much broader class of models. Al- 
though our Monte Carlo calculations were only for the 
case of Gaussian random field potentials, our analytic 
arguments above clearly hold for any generic dimension- 
less inflaton potential /, and also independently of the 
initial conditions. Specifically, in the high-energy limit 
rrih 2> fh, any messy and complicated inflaton poten- 
tial predicts the same values of (n s , r, n t , a, e, 77, £2) as the 
simple V(4>) oc <j) 2 model as long as it is not fine-tuned 
to have its second derivative vanishing in a large fraction 
of its minima that <f) can roll into, i.e., as long as generic 
minima look locally like parabolas. Of the eight cos- 
mological parameters from equation (1), we thus obtain 
sharp predictions for (O tot , w, n s , a, r, n t ). In contrast, 
the probability distributions retain finite widths for p\ 
and Q, which depend on V(<f>) and V"(<j>) at the mini- 
mum rolled down into. This means that in the rrih S> fh 
limit, the only aspect of the potential V{4>) that affects 
the cosmological parameter predictions f p (p) is the two- 
dimensional distribution of (V, V") at its minima. 



E. The low energy limit rrih <S fh 
1. Motivation 

Although the above-mentioned limit rrih S> fh can 
match observational data, it has two theoretical blem- 
ishes as illustrated in Figure 1 1 : quantum gravity correc- 
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m h /m 

FIG. 11: The inflationary model space spanned by the energy 
scales m/, and m v giving the horizontal and vertical units of the 
inflaton potential. Quantum gravity corrections may disfavor mod- 
els with > rripi PS 5m (to the right of the vertical line). The 
physically natural expectation is m v ~ m^ (green diagonal band), 
so one may hope to find a working model at the star, in the inter- 
section of this band with the cyan horizontal band that gives the 
correct fluctuation amplitude Q ~ 2 X 10 — 5 when m h ~ m P i. We 
find that such hopes are dashed by the fact that m v must be low- 
ered when mh is decreased (blue/dark grey band for Measure B), 
to offset the Q-increase caused by the preference for tiny e. One 
might naively expect quantum diffusion to be important in the yel- 
low/light grey region, but the drop in (e) as m h — > and the drop 
in (V(<^>ettd)) as m h ~ * 00 make diffusion important in the hatched 
region instead. 



tions may render mh 3> fh unnatural, and the discrep- 
ancy nth 3> m v between the two mass scales appears con- 
trived. A natural expectation would be m v ~ mh ~ m 
for some energy scale m corresponding to the physics re- 
sponsible for inflation, say the GUT scale or the string 
scale. 8 

We find that for rrih ~ fnp\, the observed fluctuation 
amplitude Q ~ 2 x 10~ 5 is reproduced if m/j ~ 0.007m, 
which suggests exploring whether inflation with = 



In the context of some supersymmetric models, the "natural" ex- 
pectation in Figure 11 is different and perhaps more compatible 
with what is observed. The reason is that in a supersymmet- 
ric theory, one naturally has two different scales: a high scale 
(say m) in the ultraviolet and the SUSY breaking scale M. In 
many supersymmetric theories, the natural scales for the poten- 
tial equation (32) of so-called moduli fields are m v = M and 
mh = m, typically with m v m/, [79]. 




100 1000 10" 
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FIG. 12: The three histograms show the distributions of the total 
number of e=foldings iVtot for high-energy (m/, — 30m, top) inter- 
mediate energy (mj, = 5m, middle) and low-energy (m/, = 0.5m, 
bottom) examples with Measure B. 



m v ~ 0.007m (five-pointed star in Figure 11) makes pre- 
dictions in agreement with observation. We will see that 
this idea, which was one of the original motivations for 
the calculations in this paper, fails. However, we will find 
that it clears one of the most challenging hurdles and 
also gives predictions largely independent of the details 
of V(4>), so it will be interesting to see whether future 
work can resuscitate it in some modified form. 



2. The flatness constraint 

Figures 9 and 10 show that as we reduce mt from the 
Planck scale to lower values, the success rate drops. The 
reason for this is obvious from equations (35) and (36): 
e and \rj\ exceed unity at most ^-values when -C fh, 
so that we obtain substantial inflation only when getting 
lucky and starting in a patch where e and \r]\ are anoma- 
lously small. The smaller mh is, the luckier we need to 
get and the lower the success rate will be. This well- 
known fact that very wiggly potentials are less likely to 
support inflation was identified early on [80-84]. 

Given that inflation is intrinsically unlikely in this low- 
energy limit, the observed flatness constraint |fi to t — 1| ^ 
0.02 [1, 2] becomes a challenging hurdle that can poten- 
tially rule out models with <C fh. If galaxy formation 
requires only N to t ^ 55 (see Section VI) and inflation 
is unlikely, is it not highly unlikely to get the additional 
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e-foldings required to satisfy the observational flatness 
constraint, A tot > 57? 

For Measure A, the answer is automatically no, since 
it automatically gives A to t = oo, so let us limit our dis- 
cussion in this subsection to Measure B. Perhaps sur- 
prisingly, Figure 10 shows that the answer is no even for 
Measure B, with low energy inflation generally passing 
this observational test. In Figure 10, the number of ad- 
ditional non-required c-foldings (-/Vbcfore) is indeed seen 
to drop from <~ 10 5 to <~ 10 1 — 10 2 as m h is reduced 
from 10 3 m down to the Planck scale, but then more or 
less stops dropping as mh is further reduced. Figure 8 



1 



-2Af b oforo 



illustrates that the histogram for |f2tot 
is almost identical for the mh = 5m w mpi case and 
that with 5 times smaller mh, typically giving values 
like ft to t = 1 ± 10~ 30 . Why is this? The basic reason 
is that, as seen in Figure 12, the A to t-distribution for 
mh <C mpi is roughly a power law /at (Atot) ~ ^tot 
a rather modest slope p <~ — 3 that steepens only slightly 
when mh is further decreased (the plotted distribution of 
lg N tot thus has slope A^ 1 ). The m,h = 0.2m case is not 
plotted in Figure 12, since it is virtually identical to the 
mh = 0.5m case shown. This means that the probabil- 
ity that the total number of e-foldings A to t exceeds some 
given value A* satisfies P(N tot < A*) oc Nl +1 ~ l/A 2 , 
so that the probability of satisfying the observed flatness 
constraint given that we are observing from a galaxy is 



/57\ p+1 
P(N tot > 57|A tot > 55) ~ ( — J 



■ 0.9. (45) 



In contrast, the flatness constraint would rule out 
mh <C m if the A to t-distribution were a much steeper 
power law or had an exponential cutoff, so let us un- 
derstand why the A tot -distribution is not steeper. If we 
rewrite the slow-roll equation (A17) as 



dN : 



m 2 (\nV)' mV27' 



(46) 



we see that there are two different ways in which we can 
get lucky and obtain large Atot when faced with tiny 
mh'. by rolling far (large A0) or by being in a region 
where the potential is extremely flat (e <C 1). As we 
will now see, our results follow from the fact that the 
latter way of getting lucky is much more likely than 
the former. To get any inflation, we need e < 1 and 
| ?7 1 < 1 at our starting point, which according to Equa- 
tions (35) and (35) requires |(ln/)'| < m h /m <C 1 and 
\f"\/f ~ (mh/m) 2 <C 1. Since / and its derivatives are 
by definition of order unity, the ^-requirement is more 
stringent and will be satisfied for only about a fraction 
(mh/rh) 2 of our simulations, showing that the success 
rate must drop at least as fast as (mh/m) 2 as mh — > 
(indeed, Figure 10 shows it dropping roughly as mh 6 )- 
This argument also predicts that almost all mh <C m 
models will have inflation terminated by \r]\ exceeding 
unity while we still have e <C 1, which is borne out 
by our simulations (for instance, out of our 1410065405 



mh = 0.5m simulations, every one of the 4260 that gave 
Atot > 55 were terminated by the 77-constraint and none 
by the e-constraint). This is why the density histories 
look so extremely flat in Figure 6 (top panel). There 
is, however, an even stronger pressure towards tiny e. 
Since the slow-roll constraints must keep holding at all 
^-values that we roll through, we have to keep getting 
lucky with e and r\ over and over again as we roll along, 
and the probability of this happening drops faster than 
any power law in the rolled distance A<fi/mh as higher 
and higher derivatives in the initial Taylor expansion be- 
come important. Figure 10 shows that empirically, by 
far the most likely way to get large numbers of e-foldings 
us therefore to roll very little, by having an anomalously 
small e and hence a tiny roll distance A(f>. Figure 10 
shows that e ~ (mh/m) 6 at the lowest energies probed, 
i.e., (InV)' <~ (mh/m) 3 , to be contrasted with the expec- 
tation (InV)' <~ (mh/m)^ 1 from equation (35). In sum- 
mary, success requires getting lucky both with e (which 
happens a fraction ~ (m^/m) 4 of the time) and with 
r\ (which happens a fraction ~ (mh/m) 2 of the time), 
which together explains the empirically observed success 
rate e ~ (m^/m) 6 . 

In conclusion, we can obtain large numbers of e- 
foldings with low-energy inflation by simply by getting 
lucky with tiny |T^"(</>)| and minuscule V'((j>) at the start- 
ing point 4>. Since zero is not a special value in the prob- 
ability distributions for V and V" (as opposed to very 
large values, say, which may be exponentially unlikely), 
this is only polynomially unlikely. 



3. Predictions for the other 7 parameters 

Above we saw that despite their low success rate, 
mh <C m models are not ruled out by their predictions for 
the curvature parameter fttot- Let us now explore their 
predictions for the remaining 7 parameters. We will first 
focus on the case of Measure A, which is arguably the 
better motivated of the two since it weights by thermal- 
ized volume. We will then turn to Measure B, and find 
that its predictions and problems are qualitatively similar 
to those of Measure A. 

Above we saw that the mh ^> m limit simplified be- 
cause the SRA was so easy to satisfy: the only properties 
of V(4>) that mattered were its minima, since the observ- 
able fluctuations where produced very near them. For 
the mh 3> m limit, there is an analogous simplification: 
the SRA is now so hard to satisfy that the only proper- 
ties of V(4>) that really matter are are the maxima, since 
the observable fluctuations are produced near them. For 
Measure A, where the thermalized pockets produced are 
infinite, we know that we began by rolling off a peak 
where — 1 < 770 < 0. Generic peaks have non- vanishing 
second derivative, i.e., fy < 0, and thus look locally like 
an upside-down parabola — for this well-know case, re- 
viewed in Appendix B.2, the distance of <j) from the peak 
grows cx a' no , i.e., exponentially with the number of c- 
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foldings. As shown in Appendix B.2, inflation around 
more general peaks is usually terminated when higher 
derivatives make rj exceed unity, and for the observable 
regime N ~ 55 gives the cosmological parameter predic- 
tions (sec Appendix B) 
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Here the ~ symbols indicate scatter of order unity in 
the relations stemming from / and its derivatives. The 
most important scatter therefore comes not from this, 
but from the strong dependence on 770 which, by the 
infinite-volume requirement, is a random variable in the 
range (—1,0). Equation (47) shows that for N ~ 55, the 
four parameters e, a, r and n* will be exponentially small 
unless r]oN < 1, i.e., unless —0.02 < 770 < 0. Since equa- 
tion (36) shows that the curvature 770 is much larger at 
generic extrema (\t]\ 3> 1), we thus expect this exponen- 
tial suppression of order 98% of the time 9 . Even for the 
remaining ~ 2% of the models, these four parameters will 
be tiny, since equation (47) gives e < N^ 2 (mh/fh) 6 < 
10 -9 for mh < 0.1m. 

Figures 5, 7 and 9 show that these analytic predictions 
are borne out by our numerical calculations: for m,h <C 
to, the parameters e, a, r and n t are all near zero and 
n s asymptotes to a broad distribution with range — 1 < 
n s < 1. As suggested by equations (49) and (51), the 
distributions for InQ and n s are seen to have the same 
shape up to a sign reversal, this common shape being 
simply that of the 770-distribution since up to additive 
factors, — InQ and n s equal 2?7o and Nrjo, respectively. 

Whereas the predictions when m/, ^> fa agreed well 
will observational data for a suitably tuned m^-value, 
the predictions when mh <C m agree quite poorly with 
the observed spectral index n s w 0.98 ± 0.02, mostly 



9 Specifically, this suggests that the ^-distribution cannot have 
a characteristic scale on this comparatively tiny sub-interval 
(—1, 0), and therefore can be well approximated by a power law. 
The In Q-distribution is the convolution of the broad and fea- 
tureless power-law distribution for Nr/o with a much narrower 
distribution of width A In Q ~ 1 reflecting when exactly the /- 
factors end inflation, and will therefore have essentially the same 
power-law distribution as Nrjo over most of the range, i.e., be 
for all practical purposes determined by the distribution of peak 
curvatures rjo alone. 



being way too red. In addition, equation (49) shows that 
generic 7y -values predict Q- values that are exponentially 
large compared to the vertical energy scale m v : for the 
extreme case 770 = 1, matching the observed value Q rts 
2 x 10~ 5 would require 
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These two potential problems are clearly linked, since 
they both correspond to 770 <C 0. We will return to 
these interesting issues in Section VI, since they require 
a careful treatment of selection effects related to the Q- 
parameter. 

Although our Monte Carlo calculations were only for 
the case of Gaussian random field potentials, our analytic 
arguments above clearly hold for any generic dimension- 
less inflaton potential /, and also independently of the 
initial conditions. Specifically, in the high-energy limit 
to/j » to, any messy and complicated inflaton poten- 
tial predicts negligible values of (r, nt, a, e, £2) as long as 
it is not fine-tuned to have its second derivative vanish- 
ing in a large fraction of the maxima with 77 > — 1, i.e., 
as long as generic maxima look locally like upside-down 
parabolas (we cover the case of unbounded potentials 
in Section V G) . Of the eight cosmological parameters 
from equation (1), we thus obtain sharp predictions for 
(fltot,w, a,r,n t ). In contrast, the probability distribu- 
tions retain finite widths for n s and Q (which are highly 
correlated since they both depend on the peak curva- 
ture 770) and for pa which depends on V{4>) at the sub- 
sequent minimum. In Section VI, we will see how the 
PA-distribution can be computed analytically from the 
distribution of the other parameters, so this means that 
in the to^ ^> to limit, the only aspect of the potential 
V{4>) that significantly affects the cosmological parame- 
ter predictions / p (p) is the distribution of curvatures 77 
at the peaks. 

The situation for Measure B is somewhat more com- 
plicated since, as seen in Figure 6, it includes cases where 
the only substantial period of inflation occurred at a 
near-inflection point rather than close to maximum, and 
cases where most of the inflation occurred after our cur- 
rent horizon scale fluctuations had been generated. Fig- 
ure 8 shows that this favors n s -values exceeding unity 
and gives a broad a-distribution in poor agreement with 
the observed constraint \a\ < 0.01 [3]. However, Mea- 
sure B is seen to share a problematic broad distributions 
for InQ. The reason is the same for both measures: 
Long-lasting inflation requires a tiny e which can give 
huge Q oc Qt/y/e, and the value of the e-parameter 55 
e-foldings before the end of inflation depends sensitively 
on the second (and for Measure B also higher) derivatives 
of the potential. 
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Scalar spectral index n s 

FIG. 13: Constraints and predictions in the (n s , r)-plane. The 
nested shaded regions are ruled out at 95% confidence from WMAP 
alone [1], when adding SDSS galaxy clustering information [2, 86] 
and when also adding SDSS Lyman a Forest information [3]. The 
green/grey points are the predictions from our simulations with 
mfc = mpi and Measure A. The two dotted curves delimit the 
three classes of inflation models known as small-field, large-field 
and hybrid models. Some single-field inflation models make highly 
specific predictions in this plane as indicated. From top to bottom, 
the figure shows the predictions for V((j>) <x 6 (line segment; ruled 
out), V(<j>) oc </> 4 (star; ruled out), V(<f>) oc <fp (line segment, 2> 
mpi; still allowed), and V(<f>) oc 1 — (</>/0*) 4 (horizontal line segment 
with r ~ 10~ 3 ; still allowed), assuming N = 64 for the c/> 4 model 
and 50 < N < 60 for the others as per [36]. 



F. The intermediate case mj, ~ m: Planck-scale 
inflation 

Above we saw that models with both mh 3> fh and 
TUh <C fh face severe problems, the former theoretically 
and the latter observationally. Let us therefore briefly 
discuss the third option, mh, ~ fh, i.e., that the horizon- 
tal scale of the inflaton potential is of order the Planck 
scale. 

This is the most complicated of the three cases: the 
predictions f p (p) will depend on the full details V(<p), not 
merely on the peaks or troughs, and none of the above- 
mentioned analytic approximations hold. The only rela- 
tion between the observables that still holds is therefore 
the familiar slow-roll consistency relationship n t = — r/8. 
To make interesting predictions for this case, it is thus 
crucial to have a physically motivated model for V(<j>), 
which we currently lack. We will therefore limit this sec- 
tion to one interesting feature which should be rather 
generic. 



Figure 13 compares the predictions the spectral index 
n s and the tensor-to-scalar ratio r for about 10 3 of our 
to/i = mpi simulations with Measure A. The observa- 
tional constraints will continue to improve dramatically 
in coming years — for instance, the completed SDSS sur- 
vey combined with the Planck CMB measurements have 
been forecast [85] to measure both of these two parame- 
ters to an accuracy of about 0.01. For n s , this means that 
all three of the cases we have considered (mh <C fh, ~ 
fh and m/, ^> fh) generically predict departures from 
scale-invariance (n s — 1) that should soon be detectable. 
As regards gravitational waves, the r-predictions are seen 
to scatter up to r-values greatly exceeding 0.01, with the 
mean being at the readily detectable level r « 0.03. This 
bodes well for upcoming CMB polarization missions aim- 
ing to detect B-polarization. 

Both of these qualitative predictions for the ~ mpi 
case (good prospects for detecting gravitational waves 
and lack of scale invariance) are rather generic, since we 
have seen that they with a greater margin (n s = 0.96 and 
r w 0.15) for essentially any potential when 3> mp\. 



G. The effect of quantum diffusion and unbounded 
potentials 

Our calculations above ignored the effect of quantum 
diffusion [32] and did not include unbounded potentials 
where V — > oo — let us briefly comment on these two 
issues. 

For these calculations, diffusion makes essentially no 
difference for Measure B and only a minor difference for 
Measure A, changing the half-basin weights by factors of 
order unity. Specifically, quantum diffusion is negligible 
on parts of the inflaton potential producing fluctuations 
Q < 1 [32] . For the calculations in this section, diffusion 
is therefore important only right near the top of peaks 
in V(<p) (where we start with Measure A), and will then 
have the effect of giving roughly equal odds for rolling 
down from this peak to the left and to the right. In other 
words, including it numerically would merely readjust the 
weights of each pair of half-basins separated by a peak 
by a factor of order unity, making the two weights equal 
rather than proportional to their length. 

Having the potential blow up (V — > oo) in places is ei- 
ther irrelevant or observationally ruled out, depending 
on the horizontal mass scale m^. If rrih 3> fh, then 
as we saw above, the observable cosmological parame- 
ters depend only the behavior of V near its minima and 
are unaffected by how the field rolled there (except in so 
far as it affects the inter-basin weighting). If <C fh, 
then 55-e-foldings before the end of inflation, <f> was most 
likely either exponentially close to a maximum of V or 
extremely high up a V — > oo hill. Although our lack of a 
theory of quantum gravity precludes us from calculating 
predictions when V((j>) ^> fh 4 , even the less extreme case 
V(4>) ~ m 4 gives Qt ~ 1, i.e., tensor modes producing 
large-scale unpolarized CMB fluctuations about 5 orders 
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or magnitude larger than the current observational lim- 
its. 



VI. CONDITIONING ON REFERENCE 
OBJECTS 

Above we computed probability distributions for the 
cosmological parameter vector p for various potentials 
and measures. However, as discussed in Section III, these 
distributions (Figure 7 and Figure 8) are not the final the- 
oretical predictions that we are entitled to confront with 
cosmological observation, since they do not include selec- 
tion effects and hence lack the appropriate weight func- 
tion w(p) in equation (29). They implicitly used protons 
or thermalizcd physical volume as reference objects, but 
the vantage points from which we make our observations 
are not a random protons but highly unusual protons (on 
planets, etc.). Instead, we can think of distributions like 
those in Figure 7 and Figure 8 as the raw probability 
distributions that come out of inflation, to be modulated 
by conditioning on whatever is appropriate to condition 
on [46, 47]. To make this more explicit, let us rewrite 
equation (29) as 



Wiv-factor reflects the fact that long-lived galaxies are 
only formed if the number of e-foldings exceeds 



/ P (p) = /; nf (pMp), 



(56) 



defining fp (p) = ((/ p (p; 3>, V))q)v as the raw parame- 
ter distribution emerging from inflation. These two parts 
of the problem (computing /™ f and computing w) thus 
decouple completely, simply being multiplied at the end. 
Likewise, the resulting distribution f p for some rudimen- 
tary reference objects (protons, say) simply gets multi- 
plied by an additional weight function if we compute the 
distribution f p for more complex reference objects made 
of the rudimentary ones (galaxies, say). We will there- 
fore explore only a limited set of reference objects in this 
section, making various simplifying approximations, leav- 
ing a more thorough exploration of weight functions for 
future work. Specifically, we will explore various combi- 
nations of the following weight factors: 



«W(P) oc e 3Wtot 
Whaio(p) ~ erfc 



n 1 V 3 



(57) 
(58) 
(59) 



w N (p) w 9[N tot - AT ga i(Q)], 
w q(p) = broad function around Q <~ 1CT 5 , (60) 

where £ = £ b + £c + & « 3.3 x 1(T 28 (sec Table 1) is 
the matter-to-photon ratio in Planck units. w\roi(p) is 
simply the factor by which inflation expands the volume. 
^haio(p) is the fraction of all protons that end up grav- 
itationally bound in nonlinear objects (in dark matter 
halos), so including the Wh a i -factor in w(p) corresponds 
to narrowing the class of reference objects to gravitation- 
ally bound protons. If the reference objects are galaxies, 
then the WN-i&ctor and the lower cutoff in the wg-factor 
are arguably required as well. As described below, the 



A^gai(Q) « In 



' 1/6 1/12 * 
Pinfl ^reheat 
£Ql/2 



(61) 



- otherwise, pre-inflationary fluctuations cause dark 
matter halos to be engulfed by black holes. The lower 
cutoff in the wq (<2)-factor reflects gas cooling physics 
related to galaxy formation and the upper cutoff incor- 
porates effects that are relevant if the reference objects 
are stable planets or observers [87]. 

Since the inflationary parameters N totl Q and p\ upon 
which equations (57)- (60) depend are correlated with our 
other five parameters, applying these selection effects 
modifies the observational predictions for all eight in- 
flationary parameters in Table 1 (and also for £ if this 
can vary from basin to basin). Arguably our most strik- 
ing result in this section will be that when taking the 
selection effects into account, broad classes of inflation 
models are ruled out by the low observed CMB fluctua- 
tion amplitude Q ~ 2 x 10~ 5 — we will refer to this as 
the smoothness problem below. 

In the next four subsections, we will derive and dis- 
cuss the weight factors in equations (57)-(60). We will 
then discuss the implications in sections VIE and VI F. 
Figure 1 shows an example of the final outcome of these 
calculations: inflationary predictions with these weight 
factors (selection effects) included. 



A. Volume weighting and conditioning on protons 

Early on, the number density of quarks (later to end 
up bound into protons and neutrons) is comparable and 
proportional to the number density of photons n 1 . Then 
matter and antimatter annihilate, leaving only a small 
residual proton density n p — r\n^, where the baryon-to- 
photon ratio r\ k> 6 x 10~ 10 (in terms of the parameters in 
Table 1, t] = ^\,/m p ). As long as the physics of baryoge- 
nesis that determines r] has nothing to do with inflation, 
the number of protons N p in a given comoving region of 
space that had volume V re h, temperature T ro h and den- 
sity p le h at the end of reheating is thus given by 

N p = V rch n p oc V rch n 7 cx V Teil T? eh cx V rc hPreh- ( 62 ) 

If the reheat density is independent of the basin into 
which we roll down (hence constant across our ensem- 
ble) and p cx a~ 3 during reheating as usual [36], then 
the number of protons is thus proportional to the phys- 
ical volume at the end of inflation. The key question is 
whether this justifies including the volume factor of equa- 
tion (57) in the weight function w{p) of equation (29). 

When comparing finite volumes, the answer appears 
to be a clear yes, since there will be only a finite number 
of protons and hence no ambiguities associate with their 
ordering in equation (26). This implies that when com- 
paring finite and infinite pockets, the finite ones get zero 
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statistical weight since the infinite ones contain almost 
all reference objects. In terms of the previous section, it 
argues for Measure A over Measure B. When comparing 
two different infinite volumes, the answer is less clear, 
and is equivalent to specifying the inter-pocket weight- 
ing. Are all infinite volumes equal, or are some more 
infinite than others? We will return to this important 
question in Section VII. 



B. Whaio: conditioning on halos 

This particular selection effect (that we are observing 
from a region of space where gravitationally bound ob- 
ject have been able to form) was the one used to obtain 
the classic anthropic upper bounds on the dark energy 
density p\ [61, 64, 67-71]. Its inclusion is clearly not op- 
tional when testing theories where Q, pa or £ vary, and 
it is more elegant than many other selection effects in 
that it is physically clean, involving only well-understood 
gravitational physics and no uncertainties related to the 
formation of galaxies, planets of observers. Any future 
definition of an observer is likely to require the inclusion 
of the Whaio-factor, since no complex objects whatsoever 
(let alone observers) ever form in the linear regime. 

The first mass scales M (if any) to go nonlinear arc 
the smallest, which enter the horizon well before matter- 
radiation equality with a fluctuation amplitude <7*(M) ~ 
Q and then grow over time due to gravitational instability 
as (see Appendix A of [88]) 



a{M) 



1 + 



r m 
PA 



1/3 



Pa 

Pm 



a.(M), (63) 



where the dimensionlcss function 



1/3 



X 

G2, 



0.7!).', 



(64) 



describes how fluctuations grow as the cosmic scale factor 
a as long as dark energy is negligible (G\{x) ~ x 1 / 3, = 

(pa/ pm) 1 ^ 3 p m oc a for and then asymptote 

to a constant value as a — > oo and dark energy dominates: 
Ga(x) — > Goo as x — > oo, where 

5r(I)r(f) 



G r , 



3V^ 



1.43728. 



(65) 



Equation (64) assumes negligibly small neutrino masses, 
PA > and a flat Universe — we will discuss the pA < 
case below and include the effect of spatial curvature in 
Section VI C 2. Neutrino masses lower the pA-predictions 
marginally [88, 89], but within the stringent current ob- 
servational limits [3], they have negligible effect on the 
qualitative results of this paper. 

In Planck units, the matter density at matter-radiation 
equality is 



Q oq _ n cq<r _ 2 C(3) y 3 



£• 



(66) 



where the Riemann zeta function £(3) w 1.20206 and the 
matter-radiation equality temperature is 



30£(3) 



21 /£\ 4/3 



£« 0.220189£, (67) 



where the second term in square brackets reflects the con- 
tributions from neutrinos to the radiation density. Equa- 
tion (63) thus shows that the net fluctuation growth is 
controlled by the factor 



t 

PA 



1/3 



— w 0.13751 



£ 4/3 

— T73 • 

Pa 



(68) 



observed to be 3215 ± 639 [2, 88], which is the factor 
by which the universe expanded between matter domi- 
nation (when growth effectively begins) and dark energy 
domination (when growth effectively ends). 

We approximate the fraction of matter collapsed into 
dark matter halos by the standard Press-Schechter for- 
malism [90], which gives 



Whaio(p) = erfc 



5 C 



V2a(M) 



erfc 



2.7046py 3 

£ 4 /%*(M) 



, (69) 



where 5 C = (9/5)2- 2 / 3 w 1.63 (not 1.69) is the fluctu- 
ation threshold corresponding to the infinite future (see 
[68] and appendix D.4 of [88]). In the second step of equa- 
tion (69), we used equations (64), (65) and (68) and took 
the limit t — > oo, i.e., p m — ► 0. To derive equation (58) 
from equation (69), all that remains is to relate the initial 
fluctuation amplitude a(M) to Q. For the small scales 
relevant to galaxy halos, the ratio <r(M)/Q is a constant 
of order unity which depends very weakly (roughly log- 
arithmically) on the mass scale M. To be specific, we 
use the observed ratio a{M)/Q w 29 corresponding to a 
galactic halo scale M = 10 12 M Q [2, 88]. This particular 
choice of scale is irrelevant to the qualitative conclusions 
we will draw in this paper, since slight variations in the 
numerical prefactor 0.1 in equation (58) are dwarfed by 
the order-of-magnitude changes in pA and Q. 

Equation (58) applies to the case pa > and gives 
Whaio ~ if pa > £ 4 Q 3 - If PA < -£ 4 Q 3 , then space will 
recollapse before any reference objects have had time to 
form [61, 64, 67-71]. Although the physics is more com- 
plicated than for the pa > case (since fluctuations do 
have time to go nonlinear before the Big Crunch), there 
is no indication that the constraints on negative values of 
Pa are much weaker than those on positive values. This 
implies that 



Whalo(p) 



1 if Ipa|«£ 4 Q 3 , 
Oif Ipa|»£ 4 Q 3 , 



(70) 



in Planck units, which is all we need for the qualitative 
conclusions of the present paper. The origin of this result 
is that the parameter combination pa/£ 4 Q 3 is simply the 
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ratio of the dark energy density to the matter density at 
the epoch when halos would form if p^ = 0. Regardless 
of the sign of p\ 1 dark energy domination is when time 
starts running out for the formation of reference objects: 
either quite literally for the p\ < case (because space 
recollapses) or for the p\ > case (because fluctuation 
growth stops). 

C. wn and wq: conditioning on galaxies 



where R = Preheat /pinfl equals unity if reheating is in- 
stantaneous and decreases if the duration of the reheat- 
ing epoch grows. For canonical values £ w 3.3 x 10 -28 , 
Q w 2 x 1(T 5 and p inR ~ (10 16 GeV) 4 - 1CT 12 , we thus 
obtain A ga i w 55 for instant reheating (R = 1). Lowering 
the reheat energy to the extreme value 1 TeV (arguably 
about the lowest energy that is still observationally al- 
lowed) gives iVgai ~ 65 and, conversely, lowering the in- 
flation energy scale reduces A ga i. 



Above we computed the weight factor Whaio that en- 
tered when the reference objects were dark matter halos. 
Let us now compute the additional weight factors that 
enter when the reference objects are galaxies. We will 
see that they depend primarily on N tot and Q. 



1. When galaxy-scale fluctuations were generated 

Consider a comoving Hubble volume (either ours or a 
smaller one) that left the horizon at a ex it during infla- 
tion (open triangle in Figure 3, say) and reentered at ciq 
during matter domination (filled triangle, say). As dis- 
cussed in Section II and Appendix A, these two points 
are connected by af ) p{ao) = o 2 xit p(a e xit). By making the 



standard assumptions that p 



exit 

OC 



-3 



during matter dom- 



ination, p 



during radiation domination, p 



oc 



during reheating and p °S a during inflation (the five 
above relations are illustrated in Figure 3 by lines of slope 
—3, —4, —3, and —2), one readily derives the follow- 
ing expression for the number of c-foldings N between 
horizon exit and the end of inflation [36] : 



AT = ln^l«ln 

f^exit 



1/6 1/12 
Pinfl Preheat / °0 
1/4 
Peq 



1/2 



(71) 



Working in Planck units, combining equation (66) with 
equation (67) shows that p cq ~ £ 4 . We use the symbol ~ 
to indicate that we are ignoring numerical factors of order 
unity for simplicity, since equation (71) shows that they 
affect No only logarithmically. The number of e-foldings 
N eq between horizon exit of the matter-radiation equality 
scale (open square in Figure 3) and the end of inflation 
(five-pointed star) is thus simply 



N 



eq 



/ 1/6 1/12 \ 
m I PinR Preheat \ 



(72) 



Fluctuations that entered the horizon before matter- 
radiation equality grew during matter domination with 
amplitude of order (a/a oq )Q, so the first nonlinear ob- 
jects formed at the epoch a ~ a cq /Q. Since galaxies 
cannot form before this time, the number of e-foldings 
Agai between horizon exit of the galaxy scale and the 
end of inflation is therefore at least 



N 



gal 



In 



, o 1/4 R 1 / 12 ^ 

Pinfl-"- 



(73) 



2. The black hole constraint on galaxy formation 

In the context of galaxy formation, there are two 
sources of fluctuations, one good and one bad: 

1. Inflation-generated fluctuations that are of order Q 
when they enter the horizon. 

2. Pre-inflationary fluctuations that we assume to be 
of order unity when they enter the horizon. 

Consider what happens in a spatial region that inflated 
by only a finite number of e-foldings N to t that exceeds 
Acq ~ 50, say. Since inflation diluted away sub-horizon 
pre-inflationary fluctuations, the first fluctuations to en- 
ter the horizon after matter-radiation equality have am- 
plitude Q, then grow oc a due to gravitational instability 
as long as pa is negligible, perhaps forming galaxies if 
they are able to grow enough to go nonlinear. If p\ = 0, 
then the horizon will keep growing until, eventually, pre- 
inflationary fluctuations of order unity start entering the 
horizon. This is bad news, since once the current horizon 
scale has gone nonlinear, the horizon volume lies within 
its own Schwarzschild radius and will recollapse into a gi- 
ant black hole on the local Hubble timescale, i.e., roughly 
by the time the cosmic age has doubled. 10 Even if the 
first order-unity fluctuations to enter the horizon happen 
to be negative and void-like, still larger-scale fluctuations 
will keep entering and ensure that our local spatial patch 
is engulfed by a black hole before long. In summary, 
horizon-sized black holes form at the epoch 



«bh ~ a eq e 



2(JV tot -JV eq ) 



(74) 



when pre-inflationary fluctuations start entering the hori- 
zon. 

If our reference objects arc merely nonlinear struc- 
tures, then fluctuations of the above-mentioned types 1 
and 2 are equally useful for producing them, and there is 
no constraint whatsoever on A to t , the duration of infla- 
tion. Indeed, no inflation at all is needed, since objects 



10 Note that fluctuations of order unity in a flat Friedman- 
Robertson- Walker background only form black holes if they are 
on or above the current horizon scale. In contrast, fluctuations 
with Q <K 1 grow to order unity only much later, when the cur- 
rent horizon scale is vastly larger than the fluctuation scale. 
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would be nonlinear from the outset. In contrast, refer- 
ence objects such as galaxies, planets and observers all 
require extra time to form. Specifically, they require a 
period of peace and quiet after their parent dark mat- 
ter halo has gone nonlinear, and will not have time to 
form if their parent halo is promptly swallowed by a 
black hole. Such reference objects will therefore only 
form if halos form substantially before the epoch a^, 
i.e., if N tot > N ga \(Q). This constraint on the duration 
of inflation is embodied by equation (59). 

None of the qualitative conclusions in this paper de- 
pend on whether -/V ga i = 50, 60, or some other number of 
order 10 2 . Since the dependence of A ga i on Q is only log- 
arithmic, i.e., much weaker than that of the Whaio factor 
of equation (58), we therefore made the crude approxi- 
mation iVg a i = 55 for simplicity in Section V. 



3. wq : The cooling constraint on galaxy formation 

To end up in a galaxy, a proton must first fall into a 
dark matter halo (a fraction tUhaio factor does this), then 
avoid promptly being engulfed by a black whole (a frac- 
tion wn avoids this). Next, the gas to which it belongs 
needs to cool efficiently enough to be able to contract, be- 
come self-gravitating and form stars (see, e.g., [91-94]). 
This problem is analyzed in terms of fundamental physi- 
cal constants in [87]. If we decrease Q, then dark matter 
halos form later and have lower escape velocities, hence 
heating infalling gas to lower temperatures where cooling 
is less efficient, causing the weight function wq(Q) to fall 
off towards very low Q- values. After changing variables 
to reflect the notation in this paper, equation (11) in [87] 
states that wq(Q) ~ unless 

Q > 0.2a- 1 ln[ a - 2 ]- 16 / 9 m4/3 ?? -2/3 r 2/3 _ 1Q - 6j (?5) 

where a w 1/137.036 is the fine structure constant and 
Top ss 4.2 x 10~ 23 is the electron mass in Planck units. 



D. wq: conditioning on stable planets 

If the reference objects are observers who, like us, have 
required billions of years of evolution on a planet, we need 
to include an additional weight factor wq(Q) incorporat- 
ing the fraction of the protons in galaxies that end up in 
planets in long-lived stable orbits. As discussed in [87], 
if we increase Q, then galaxies form earlier with higher 
stellar density, increasing the rate of catastrophic orbit 
disruptions from near encounters with other stars. Equa- 
tion (19) in [87] states that wq(Q) « unless 



Q < 10" 



(76) 



In summary, combining equations (75) and (76) suggests 
that wq(Q) ~ unless 10~ 4 < Q < 10~ 6 . However, as 
detailed in [87], the physical assumptions underlying the 
calculations of wq are not nearly as clean and clear-cut as 



those underlying, say, the Whakrfactor, so this conclusion 
should be taken with a healthy grain of salt, and deserves 
further investigation. 



E. Predictions for the horizon problem 

In the four preceding subsections, we derived the 
weight factors in equations (57)- (60). Let us now explore 
their implications, starting in this section with equa- 
tion (59). 

Inflation solves the horizon problem if A to t ^ Ao, so 
that all fluctuations that have entered our horizon so far 
have inflationary rather than pre-inflationary origin. In 
Section A 2, we saw that if inflation lasts for an additional 
few e-foldings, it solves the flatness problem as well: since 
^before = A tot - N , we found that A be fore Z solves the 
horizon problem and Abefore <; 3 reproduces the observed 
spatial flatness. 

One of the original motivations for the early work on 
inflation was that it could solve the horizon problem and 
the flatness problem [80-84] (for certain potentials). To 
what extent can we make the prediction that inflation 
does solve these problems? 

When the total number of c-foldings A to t varies across 
the ensemble, it will generally exceed Ao in some spatial 
regions but not in others. For observations made from 
galaxies, we automatically have -/V > A ga i, but the ex- 
istence of galaxies guarantees only that A tc ,t > A ga i , not 
that A to t also exceeds the larger quantity A correspond- 
ing to the current observation time. The probability of 
observing the horizon problem to be solved is thus the 
probability of observing A > A given that TV > A ga i, 
i.e., 



Phorizon — 

P(A tot > ^0|Atot > A gal ) = 



P(A tot > No) 



P(N tot > A gal )' 
(77) 

In Section V, we found that Measure A predicted A tot = 
oo, automatically solving the horizon and flatness prob- 
lems, whereas our other measure gave a probability dis- 
tribution for A falling off roughly as A t ~7 with 7 m 3 in 
the worst-case scenario rrih <C ffi (plotted in Figure 12). 
Equation (71) shows that A — A ga i = i ln(l + z ga i) ~ 1, 
where z ga i is the galaxy formation rcdshift. For A ga i = 
55, equation (77) therefore gives 



Pho 



1 with volume-weighting, 



No 



~ 1 - i '^r ~ 95% otherwise- 

(78) 

For either measure, inflation thus predicts that the 
horizon problem is solved from the vantage point of 
most galaxies. If we were to retrodict whether our 
own horizon-scale fluctuations were inflationary rather 
than pre-inflationary, varying z ga i over the observation- 
ally rather extreme range 1 < z ga i < 100 would give 
87% < Phorizon < 98%, so our qualitative conclusion that 
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^horizon ~ 1 is independent of uncertainties regarding 
what to take as the present epoch. We can also make this 
argument in a less anthropocentric way, noting that since 
the Whaio-factor tends to postpone dark energy domina- 
tion tends to only slightly after galaxy formation, the 
horizon stops growing and precludes observers from ob- 
serving more than a few e-foldings beyond N ga \ no matter 
how long they wait. 

In conclusion, none of the inflation models we have 
considered are ruled out by the horizon problem. In Sec- 
tion V, we saw that they were not ruled out by the flat- 
ness problem either. 



F. The smoothness problem and other implications 

If a physics problem involves two vastly different scales, 
this fact can often be used to simplify it. We will now see 
that precisely this happens in our inflationary prediction 
problem, because the inflation density vastly exceeds the 
dark energy density £ 4 Q 3 m P1 allowed by equation (70). 
Conveniently, this will allow us to deal analytically with 
the pA-parameter and dispense with it once and for all. 

In Section V, we saw that inflation typically ended at a 
density p on( j ~ m„ 4 for rrih ^ fh and p cn( j <~ (m^/m) m v 
for nih ?J fh. For any relevant values of Q, and m v , 
both of these two densities vastly exceed the density scale 
£ 4 Q 3 TOp[ <~ lO -110 ^ 3 ?™ 4 ,! above which equation (58) cuts 
off the statistical weight (no halos form) . The lower limit 
on negative pA-values from requiring reference objects to 
form before space recollapses is of comparable magnitude. 
This means that the only part of the probability distribu- 
tion for p A = y(0 st op) = m^f^stop/rrih) that matters 
is that in a region around zero which is way smaller in 
magnitude than any relevant physical scales. Since the 
dimensionless inflaton potential / by definition varies of 
order unity and zero is not in any way a special value for 
a random minimum of a complicated function, it follows 
that we can for all practical purposes take pa to have a 
uniform distribution, uncorrelated with all the other pa- 
rameters. This prediction is confirmed by our numerical 
results and illustrated by the pA-histograms in figures (7) 
and (8). 



1. Marginalizing over pA 

Since Whaio is the only one of the weight functions from 
equations (57)-(60) that depends on pa, this result means 
that we can once and for all marginalize over p\ analyt- 
ically and eliminate it from the discussion of our other 
cosmological parameters. The resulting parameter dis- 



tribution observed from a random halo is 

/oo 
f P nf (p)w halo (p)dp A 
-oo 

« /^4 nf (P)|,A=0ff(^)^A 

cx £ 4 Q 3 4 nf (p)| PA =o. (79) 

On the second line, we used the above uniform distri- 
bution result to approximate the first factor /p" f (p) by 
its value at pa = because it is pA-indcpcndent when- 
ever the second factor is appreciably nonzero. 11 Here 
we also reexpressed Whaio(p) as <7(pa/£ 4 Q 3 ) to empha- 
size that the final result in equation (79) applies to much 
more general weight functions u> halo (p) than of that of 
equation (58) — which corresponds to the particular case 
g{x) = erfc (O.lx 1 / 3 ). Replacing the Press-Schechter for- 
mula by the more accurate approximation of Sheth & 
Tormen [95] would merely modify the functional form of 
g(x), not the final result in equation (79). More impor- 
tantly, modifying Whaio(p) to explicitly give zero weight 
to models where space recollapses before reference ob- 
jects form (which is necessary to prevent the integral in 
equation (79) from diverging as pa — > — oo) would again 
simply modify the functional form of g(x), this time to 
ensure that the condition g(x) w for |a;| » 1 held also 
for negative x, as per equation (70). 

The Q 3 -factor on the last line of equation (79) poses 
great difficulties for many inflation models, since it can 
easily overpower a slight inflationary preference for low 
Q- values in / lnf (p) (Figure 14) and give predictions far 
exceeding the observed value Q^2xl0 -5 . As a con- 
crete illustration of this smoothness problem, consider the 
models with nth < fh from Section V. These were the 
only models that were natural, in the sense that quantum 
gravity corrections may preclude larger m^-values. We 
found that all these models predicted a Q-distribution 
with a high tail falling off roughly as Jq{Q) oc Q _1 , cor- 
responding to a near- uniform distribution for InQ, as 
confirmed numerically in Figure 14. The basic reason for 
this was that Q <~ V z / 2 /\V'\, and that V was often ex- 
ponentially small. Equation (79) now shows that these 
models predict a rising tail, with a random halo seeing a 
distribution /q(p) oc Q 2 for large Q. Quantum diffusion 



Because our simulations predict the pA-distribution to be uni- 
form over a much broader range (|pa| "C Pend) than the range 
forming halos, it would be numerically inconvenient to apply 
equation (58) to our simulations: it would give a success rate 
<C 10 -100 and require many Hubble times worth of computer 
simulating to obtain even a single successful model. Requiring 
PA = as in equation (79) is of course even less convenient 
numerically. In practice, we therefore use the result of equa- 
tion (79) with the constraint pa = replaced by < pa < p\, 
where p A <C p en d- We choose the cutoff p A sufficiently small that 
reducing it further has no effect on our predicted distributions, 
but large enough that we retain a statistically large sample of 
models. 
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FIG. 14: The smoothness problem is seen to afflict quantum- 
gravitationally "natural" inflation models with m^ <C mpj Ri 5m, 
predicting Q way above the observed value Q ~ 2 X 10 — 5 since 
PA-marginalization multiplies the above distributions by an extra 
Q 3 -factor. The three histograms are for Measure A and show the 
distribution of lg Q for pa ~ 0, marginalized over all other cos- 
mological parameters. The = 0.5m distribution is seen to ap- 
proach the analytic m h — * prediction of a uniform distribution 
in the high lgQ tail, and the ra/, = 30m curve is seen to be ap- 
proximated by the analytic m/, — > oo prediction <x Q 2 e _ ^/ < ^*) 
for Q» ~ 5m v 2 /m h near the peak (dashed line). 



cuts off this distribution at Q ~ 1, so Q observed from 
a random halo would be of order unity. Equation (55) 
shows that this result is independent of the choice of the 
vertical energy scale m v unless it is down in the TeV 
range — for that extreme case, equations (49) and (51)) 
show that ijo « — 1 and n s w —1. 

All these <C fh models are thus firmly ruled out by 
the observed measurement Q~2xl0~ 5 unless they are 
rescued by the weight function wq(Q). 12 However, as 
discussed in [87] and Section VI C 3, the anthropic upper 
limits on Q from planetary orbit disruptions and other 
effects appear rather weak, both in terms of physical un- 
certainties and in terms of the actual numbers. In partic- 
ular, there is currently no compelling argument for why 
there would be much fewer observers if Q were say 3 times 
larger than observed, yet these inflation models predict 



FIG. 15: Predictions for the dark energy density pa for the toy 
model where reference objects get produced if 10 -6 < Q < 10~ 4 
and inflation predicts a power law /q(Q) oc Q n across this range. 
Peaking from left to right, the solid curves correspond to n = — 8, 
— 7, —6, —5, —4, —3, —2, —1 (heavy), 0, 1, 2 and 3, respectively. 
The short-dashed curve shows the case where inflation predicts 
Q = 2 X 10 -5 without any scatter, and the long-dashed curve 
is for the actual Q-distribution for m h = 10m, m v = 0.01m from 
Figure 1. The shaded band shows the WMAP+SDSS measurement 
from Table 1. 

the (In Q)-distribution a Q 3 so that this larger In Q- value 
would be a priori about 30 times more likely. Having In Q 
up by an order of magnitude (bumping up against the dis- 
ruption limit of [87]) would be a priori about 1000 times 
more likely. It would be much more satisfying if this 
smoothness problem could be solved by a well-motivated 
inflation model firmly predicting smaller Q-values. 

Analogously, one can argue that equation (79) rules out 
any particle physics model for the matter density param- 
eter parameter £ if it predicts an ensemble with a high-£ 
tail shallower than about £ -5 , unless galaxy formation 
can be shown to be strongly inhibited by small increases 
in c; — for a discussion of such dark matter issues, see 
[46, 96]. 



2. Predictions for pa 



The uijv-factor of equation (59) provides essentially no help — 
although it does depend on Q for one of the two measures we 
studied, the dependence is very weak because of the logarithm 
and the broad intrinsic width of the TV-distribution. 



Above we marginalized over p\ to obtain predictions 
for the other parameters. Let us now do the opposite, ob- 
taining inflationary predictions for p^. The inflationary 
prediction for the dark energy probability distribution 
/a(Pa) is obtained by marginalizing f p (p) from equa- 
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tion (56) over the other seven parameters in the p-vector, 
i.e., 

/a(p a ) « / / inf (pMp)%a - P ' A )dV (80) 



mm _ £4/0 
rA — S ^mm 



Above we found that / mf (p) was pA-independent over 
the narrow range of p-values where w(p) was nonzero. 
The only pA-dependence in equation (80) therefore comes 
from the weight function w(p), where pa enters in the 
halo term u>h a io(p) of equation (58) and gets coupled with 
Q. We can thus integrate trivially over all parameters 
except Q in equation (80), obtaining 



/a (pa) 



f (Q)wq(Q) erfc 



n 1 V3 



dQ, (81) 



where f lnl (Q) is the raw inflation prediction for Q for 
Pa = (i.e., f (p)| PA= o marginalized over all param- 
eters except Q), as plotted in Figure 14 and in the em- 
panels of figures (8) and (7). Since the values Qi in our 
simulations with pa ~ are drawn from the distribution 
f lni (Q), we can evaluate this expression in practice as 



/a (pa) ?S^2wq(Q 2 ) erfc 



(82) 



where the sum is over those simulations for which pa ~ 0. 
If inflation were to predict that Q took a specific value 
Q*-value (i.e., if f lnl (Q) = S(Q — Q*), as for a one- 
dimensional inflaton potential with only a single mini- 
mum that either was symmetric or could only be slow- 
rolled to from one side, then equation (81) would give the 
classic prediction p A < lO^Q 3 [61, 64, 67-71, 87] (more 
specifically, the dashed curve in Figure 15). However, as 
stressed in [97], degeneracies with other parameters can 
invalidate such arguments. Sure enough, in the cases we 
have considered in this paper, the distribution f inl (Q) 
has a nonzero width, and equation (81) shows that a tail 
towards large Q-values will give more statistical weight 
to higher pA-values. This important coupling between 
PA and Q was ignored in the early work predicting pa- 
distributions. It was emphasized in [87] and recently ex- 
plored numerically in [98] for some instructive toy mod- 
els, with the conclusion that it can invalidate anthropic 
explanations for pa for a broad class of primordial Q- 
distributions. [99] explores how affects this link. 

Using our results from Section V, we can for the first 
time explore this issue with Q-distributions from actual 
inflation calculations. 

Before doing this, however, let us first build some in- 
tuition for equation (81). Consider the toy model where 
/ (p) Q n f° r some power law index n and the weight 
function wq(Q) = 1 if Q m i n < Q < Qmax, vanishing 
otherwise. Limiting our attention to the case p\ > 
and performing the integral in equation (81) numeri- 
cally now gives the curves shown in Figure 15, which are 
easy to understand intuitively if we define the quantities 



and PA ax = £ 4 Qmax- Since erfc a: w 1 for 
i« 1 and erfcx w for x ^> 1, equation (81) shows that 
/a (pa) is constant if p\ <C p™ ln , vanishes if pa S> p™ ax 
and is otherwise approximated by 



/a(pa) 



c 4/3 



Q n dQ 



oc 



PA 

-max 

Pa 



(83) 



Figure 15 shows the distribution of lgpA rather than 
PA, which introduces an extra pA-factor since dp a = 
PA d In pa, making the exponent of the second term of 
equation (83) equal to (n + 4)/3. Whether n is smaller 
or larger than —4 therefore determines the qualitative 
behavior of the resulting distributions: Figure 15 shows 
that for n = 4, all values of lg pa between p™ 111 and p™ ax 
are roughly equally likely, whereas smaller or larger n 
cause the distributions to be peaked at p™" 1 and PA ax , 
respectively. There are thus only three classes of theories 
that are not ruled out by the observed pA-value (shaded 
band) : 



1. Theories where n 



-A. 



2. Theories where inflation predicts a very narrow dis- 
tribution f lni (Q) centered around or slightly below 
the observed Q-value. 

3. Theories where the weight function w(Q) is suffi- 
ciently sharply peaked around the observed Q- value 
to overpower the pull towards p™ ln or p™ ax - 

None of these classes are problem free. Class 1 involves 
inflationary fine-tuning in the sense that it requires the 
sharp inflationary prediction for lg Q to coincidentally fall 
within the range around —5 where wq(Q) is substantial. 
We argued in Section VI C 3 that wq (Q) may be rather 
broad, suggesting that models in class 3 will only work 
if n is not too far from —4. We saw in Section V that 
inflation models with a "natural" horizontal energy scale 
rah m do not fall into class 1, but rather that their raw 
predicted Q- values display a high Q tail with n w — 1 
(heavy curve in Figure 15). Such models thus tend to 
overpredict pa by orders of magnitude, which is yet an- 
other manifestation of the smoothness problem discussed 
in Section VI F. 

Figure 16 is an attempt to illustrate the smoothness 
problem and this intimate relation between pa and Q. In 
summary, we have found no compelling solution to the 
smoothness problem. High energy inflation models with 
nih 3> mpi can have nih fine-tuned to fall into class 1, 
but may flounder on quantum gravity corrections. Low 
energy inflation models with nih ^ mpi, on the other 
hand, predict (pa,Q) near the five-pointed star in Fig- 
ure 16, quite far from the observed value. If a future 
candidate inflation model provides a compelling solution 
to the smoothness problem, this will thus be a non-trivial 
and noteworthy achievement. 
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FIG. 16: The weight function w(p) (the number of reference 
objects produced per unit thermalized volume) causes important 
modulations of the raw inflationary predictions in the (pA;Q)- 
plane, with u>halo(p) gi vm g essentially no weight to the region 
marked "no halos" and w(Q) strongly downweighting regions more 
than about an order of magnitude above or below Q = 10 -5 . 
The fact that inflation typically predicts p\ to have a uniform 
distribution in this range, uncorrelated with Q, gives probability 
<x pa hi this plane, so nearly all probability lies within the diago- 
nal green/grey favored band near the halo limit. The logarithmic 
slope of the inflationary Q-distribution determines the most likely 
outcome, with din /din Q -C —4 giving the open square and 
dln/^ nf /dlnQ > -4 giving the star. The filled square shows the 
measurements from Table 1. 



G. Calculation summary 

Let us summarize the recipe that we have derived 
for computing the cosmological parameter distribution 

1. Use one of the two Monte Carlo approaches of Sec- 
tion V (depending on whether you wish to exclude 
finite pockets or not) to generate a large number of 
parameter vectors Pi drawn from the raw inflation- 
ary distribution /p" f (p), keeping only those with 

JVtot > ^gal « 55. 



2. To approximate the p\ — > limit of equation (79), 
discard all vectors with pa < or pA > p\, reduc- 
ing p\ until the results stop changing. 

3. Compute /a(pa) using equation (82). 

4. Compute the distributions for all other parameters 
as f p (p) cx w Q (Q)Q 3 f^ ni (p)i PA=a by weighting the 



points from Step 2 by WQ(Qi)Q 3 and then making 
histograms. 



In practice, Step 4 can be done either by unequal point 
weights WQ(Qi)Q^ in the histograms or by resampling, 
drawing a large number of points from Step 2 at random 
with probability cx WQ(Qi)Q 3 . We use the latter ap- 
proach. Figure 1 shows a complete worked example for 
the case of Measure A, m/, = lOmpi and m v = 0.002mpi. 



VII. ANALYTIC RESULTS 

In the two preceding sections, we computed cosmo- 
logical parameter predictions for specific potentials and 
measures, finding that unless m/, ~ Wpi, the dependence 
on the inflaton potential could be understood analyti- 
cally in terms of its peak and trough statistics. We did 
this with pocket-based orderings, so to better understand 
how the choice of measure affects the predictions, let us 
now revisit this computation using the other class of or- 
derings from Section IV: global orderings. We will find 
that simple analytic calculations allow us to rule out a 
broad class of global orderings, notably many that give 
attractor behavior and make the predictions independent 
of pre-inflationary initial conditions. 

In Section IV, we saw that a global ordering is uniquely 
defined by specifying a "time" variable to order the refer- 
ence objects by. We found that all observable quantities 
{e.g., p, T or H) are inadequate as such global time vari- 
ables, leading instead to pocket-based orderings, so let us 
now explore time variables that are not observable (ex- 
cept in differences or ratios, like t and a). 

We will focus most of our discussion on the case of or- 
dering by time t 13 , then discuss how these conclusions can 
be generalized to other global orderings. This natural- 
sounding choice is one of the first to have been explored 
in detail, and was used in some of the key early work on 
eternal inflation [8, 53, 62, 103-106]. 



13 Specifically, we define t as the time variable in synchronous co- 
ordinates where the metric is 



ds 2 = dt 2 - a(x, t) 2 d-x. 2 



(84) 



The lines of constant x in this metric are timelike geodesies cor- 
responding to the worldlines of Gedanken co-moving observers, 
and t is the proper time of these observers. This coordinate 
system remains well-defined until these geodesies start to cross, 
which happens only long after thermalization when fluctuations 
go non-linear — the coolness problem discussed below becomes 
severe long before then, so it cannot be circumvented by a cre- 
ative choice of gauge at late times. When we discuss t for a 
particle in the present epoch, the rigorously inclined reader can 
simply take this to mean its proper time, since this provides a 
well-defined ordering even after geodesic crossing. 
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A. The origin of attractor behavior 

One of the most appealing features of i-ordering is 
that it generically gives attractor behavior, allowing a 
unique parameter probability distribution / p (p) to be 
computed independently of the initial conditions. This 
is most elegantly seen using the Fokker-Planck equation 
formalism, giving the time-evolution of (</>) , the phys- 
ical volume associated with different </>-values. Thor- 
ough and detailed treatments of this are given in, e.g., 
[8, 53, 62, 103-106], so we will only summarize the results 
here and attempt to provide physical intuition for them. 
The Fokker-Planck equation takes the form 



(85) 



where L is a linear operator involving , I s * and 2 nd 
derivatives with respect to the inflaton field <fi to incor- 
porate the effects of expansion, slow-roll and quantum 
diffusion, respectively. The total volume at a given time 
is thus 



V(t) 



J v^4>,t)d d 



(86) 



which is conveniently partitioned as the sum of two con- 
tributions Vi n fl(£) and Vtherm(i) corresponding to inte- 
grating over the parts of </>-space where space is inflating 
and has thermalized, respectively. (The standard pre- 
scription is to reclassify a volume of space from inflating 
to thermalized as soon as its </>-value leaves the part of 
0-space where the slow-roll approximation is valid and to 
not evolve it further [8, 53, 62, 103-106], so that Vtherm(t) 
is simply the volume of the U-shaped hypersurfaces in 
Figure 4 that lie in the part of spacetime with time co- 
ordinate < t.) 

If the total volume V(t) is finite at some initial time 
t, then both Vinfl(i) and Vthcrm(i) will remain finite for 
all t since the expansion rate H cannot be infinite. In- 
flation is usually said to be eternal 14 if Vi n a(t) — ► oo 
as t — > oo. If the number of reference objects is propor- 
tional to the thermalized volume produced, i.e., Vtherm(i) 
(as is the case for points, protons, planets and all other 
reference objects we have discussed), then inflation pro- 
duces an infinite number of reference objects provided 
that Vtherm(t) — > oo as t — > co . Such eternal inflation 
with infinitely many reference objects occurs quite gener- 
ically, whenever the inflaton potential V(<p) contains an 
"eternal inflation region" . It is well known that there are 
three types of eternal inflation regions [48-50, 108]: 

1. A local maximum in V(<f>) where the slow- roll ap- 
proximation is valid [62, 109]. 



2. A 0-range where V((f>) is large enough for quantum 
diffusion to dominate the dynamics [32]. 

3. A local minimum in V((j>) from which <p can even- 
tually tunnel out ( "old inflation" ) . 

In the first case, the volume rolling off the peak in various 
directions is more than replenished by the expansion of 
the volume that has not yet rolled off. In the second case, 
quantum diffusion up the hill causes faster expansion that 
more than compensates for the fact that it would have 
been more likely to roll downward — this happens in <p- 
regions where the slow-roll approximation predicts Q > 
1. 

If there are multiple eternal inflation regions, the 
fastest expanding one (the one where V(d>) is largest) 
always wins, completely dominating the volume distri- 
bution V^t) as t -> oo. [8, 53, 62, 103-106]. The reason 
is that if two volumes are eternally expanding with rates 
Hi and H 2 , and Hi > H 2 , then the lower second peak 
will asymptotically contribute a volume fraction 



V 2 (t) _ V 2 {0)e 3H ^ _ V 2 (0) 
Vi{t) ~ Vi(0)e 3H ^ ~ Vi(0) 



e -3( Hl -H 2 )t ^ (8?) 



as t — > oo whatever the initial volume ratio V 2 (0)/Vi(0) 
was. 15 

Moreover, the fastest expanding eternal inflation re- 
gion of all (around </>„, say) will completely dominate 
the asymptotic t — > oo volume distribution regardless 
of the initial conditions, i.e., for any initial distribution 

(</>). This is because quantum diffusion will populate 
all 0-values and even an exponentially tiny contribution 
Vtj, (</>*) created, say, by exponentially suppressed uphill 
diffusion/tunneling at low density, will eventually over- 
power all other slower expanding regions as t — + oo just 
as described above. As a result, the asymptotic distribu- 
tion V,p((j))/V(t) will be determined simply by how long 
time it takes to roll/diffuse down from <j> t to 4>. Since 
we are evaluating the distribution at fixed time, faster is 
better. If a volume of space has the inflaton descending 
from along a trajectory (f>(t) during the time inter- 
val t\ < t < t 2 , then it will expand by H((j)(t))dt 
e- foldings while the eternal region expanded more, by 
(i 2 -*!)#(</>*) e-foldings (H(<j>) = m" 1 y/V (</>)/ '3), so the 
asymptotic probability distribution is [8, 53, 62, 103-106] 



V^{4>) = V</,((j)*) exp 



3 / [H(<f>.) - H(4>(t))]dt 



tl 



(88) 

Since this asymptotic distribution satisfies V<j,((f)) = 
3H*V<t,((fi), it can be computed directly from the Fokker- 
Planck equation (85) as the solution to the eigenvalue 



Throughout this paper, we use eternal to mean eternal to the 
future. For recent discussions of whether inflation can be eternal 
to the past, see [48-50, 100, 101, 107] and references therein. For 
an up-to-date discussion of initial conditions and whether they 
are needed, see [102]. 



Measure ambiguities are thus linked to the fact that we cannot 
reverse the order of the two limits n — > oo in equation (26) and 
t — > oo. 
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problem LV^ — SH^Vcf, that has the largest eigenvalue 
3-ff*. (If this largest eigenvalue is negative, inflation is 
not eternal.) This eigenvalue problem can be transformed 
into a manifestly Hermitean form, and has been shown 
to reduce to a Schrodinger equation [57] . 

What we care about empirically is not the probabil- 
ity distribution for </>, but the parameter probability dis- 
tribution ,fp(p). The latter is readily extracted from 
the former by keeping track of how inflation ends, i.e., 
by how the thermalized volume Vthorm(i) is distributed 
among the different basins of attraction and, for each 
basin, what the probability distribution is for the direc- 
tion from which it rolled down towards the minimum. 
A mathematically equivalent alternative to solving the 
Fokkcr-Planck equation (see below) is to evolve the in- 
flaton <p using a stochastic ordinary differential equation 
that incorporates the effects of quantum diffusion as a 
random noise term. 



B. What t-ordering predicts 

A strength of the t-foliation solution to the ordering 
problem is that it makes a specific and well-defined pre- 
diction for the parameter probability distribution ,f p {p) 
which, because of the attractor dynamics, is indepen- 
dent on initial conditions and unknown pre-infiationary 
physics. Let us now discuss what this prediction is. 

We saw above that in the t — > oo limit, all reference 
objects except for a set of measure zero are in comoving 
regions of space where <j> came down from the globally 
highest peak in the potential, V{4> jf ). There is ambiguity 
as to what happens if more than one peak has V rising to 
above the Planck scale, because we lack a theory of quan- 
tum gravity, but we will not worry about this issue here 
in light of the more pressing problems described below. 
Note that the attractor behavior makes not only the ini- 
tial conditions irrelevant, but also most of the properties 
of the infiaton potential V(<f>) — all that matters is the 
shape of the highest peak and its immediate surround- 
ings. 

Because of equation (88), the fastest way to come down 
from this peak is the most favored, and it has been 
shown [56] that the exponential preference for speed is 
so strong that the dominant descent mechanism can in 
some cases be quantum diffusion/tunneling rather than 
slow roll, giving the strange appearance that we are liv- 
ing in the center of a spherical void (at least for Planck- 
scale inflation). This preference for rapid descent can 
cause <f> to be quite near the top of the peak 55 e- 
foldings before inflation, so that the tensor amplitude 
Qt ~ V(4> r ) 1 ^ 2 /m 2 . This suggests that a generic infia- 
ton potential with many random-looking peaks will pre- 
dict Qt 3> 1CP 5 and hence be ruled out by observation 
unless V(4>) < ICr 10 m 4 - 10 9 GeV 4 for all 4>. For in- 
stance, the Gaussian random field potentials we explored 
in Section V are guaranteed to contain arbitrarily high 
peaks (exponentially rarely, but this is irrelevant), and 



are hence all ruled out if i-ordering is used. It should 
be emphasized, however, that there are infinitely many 
potentials that do not overpredict the fluctuation level 
with time ordering, most obviously popular single-basin 
ones like V(<p) = m 2 (j) 2 with m«5x 10~ 6 m. 



C. The coolness problem 

We will now present an argument for why the t- 
foliation solution to the ordering problem is observation- 
ally ruled out. This is closely related to arguments given 
in [48-50, 56] as described below. 

Consider any reference objects the number of which 
is proportional to physical thermalized volume, say pro- 
tons, planets or people. Let us compute the probability 
distribution /t(T) for the CMB temperature T observed 
at such reference objects at a given time t. We will find 
that our universe is much cooler than predicted. 

After a given volume thermalizes at time tthorm, 
it keeps expanding and its photon temperature Tcmb 
drops. We can therefore write Tcmb = f{t — ithcrm) 
for a known and monotonically decreasing function /. In 
Planck units (mpi = 1 ; used throughout this section) , we 
have 



/(At) 



At- 1 / 2 



during radiation domination, 



£ V 3 Ai 2/3 during ma tter domination, 



where £ w 3.3 x 10~ 28 is the matter-to-photon ratio from 
Table 1. (Tcmb a^ 1 at all times.) Since the number of 
refence objects is proportional to the thermalized volume, 
the fraction of them in domains thermalized by time t is 



P(t t herm < t) OC Vthcrm(*) »C 6 



3H,t 



(89) 



The fraction of the reference objects where Tcmb is below 
some given value T is therefore 

P[T CMB < T] = P[f(t - t thcim ) < T] 

= P[*therm < t - (T)\ 
ZH,[t-r\T)\ 



oc e 



oc e 



(90) 



where the inverse function At = f 1 (T) is approximately 



during radiation domination, 



£ V 2 T 3 / 2 during matter domination, 



Since the expression e _3H *^ in equation (90) is a 
perfectly regular and well-behaved function, approaching 
zero as T — > (t — > oo) and approaching unity as T — > 
Tthorm, eternal inflation with t-ordering thus makes the 
firm testable prediction that the temperature probability 
distribution is given by P[T CM b < T] = e^tf,/- 1 ^ 
During matter domination, we thus have the prediction 



P{Tcmb < T) w e 



-(£) 



-3/2 



(91) 
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where T* = ^ 1/3 (3i?)* /3 w 10 9 i?» /3 . For definiteness, 
let us take the eternal inflation expansion rate to be a 
typical inflationary rate H t <~ 10~ 5 /10~ 38 s, which gives 
T* ~ 10 38 K — the conclusions we will draw are com- 
pletely unaffected by varying H* within the observation- 
ally allowed range. This gives the spectacularly small 
probability 

P(T CMB <3K)«1(T 1056 , (92) 

and the probability of finding Tomb as low as ob- 
served by COBE/FIRAS, i.e., T C mb = 2.728 ± 0.004, 
is of course even smaller. Even if we take as a given 
that Tomb < 4K, this probability remains minuscule: 
P{T C mb < 3K|T CM B < 4K) w 10" 1055 ' 6 . Indeed, the 
probability is small even of having the CMB temperature 
a mere 10 _17 K warmer, corresponding to a one second 
younger (relative to the thermalization time) universe: 

about e~(! SCC ) H * rj 10 _1 ° 38 . Although it is easy to en- 
vision an anthropic weight factor w(p) that may slightly 
skew this distribution towards lower temperatures (by 
requiring sufficient time for evolution, say), such correc- 
tions will be clearly be way too week to overpower this 
double exponential preference for high temperature - 
indeed, before CMB temperature was measured in 1965, 
scientists saw no fundamental problem with it being say 
5K as predicted by Alpher & Herman [110]. The FIRAS 
measurement of the CMB temperature therefore rules out 
eternal inflation with time ordering at a significance level 
of about 99. 999. ..9%, where the dots denote 10 56 nines. 
Few scientific theories have had the honor of being ruled 
out with greater statistical significance. 

The root of the failure is clearly that this ordering 
rewards rapid expansion, so that at a fixed time, the 
vast majority of all reference objects formed only very 
recently. Although we focused on Tomb above to be spe- 
cific, the coolness problem clearly afflicts many other ob- 
servable parameters as well. Our universe is not only 
cooler than predicted, but also older and less dense than 
t-ordering predicts. These paradoxical implications have 
been discussed in the literature, although in contexts 
leading to somewhat less negative conclusions. This is- 
sue is at the core of the surprising findings of Linde & 
Mezhluminan [56] that we might expect to find ourselves 
in the center of a rare spherical void, and the authors 
mention this as a perhaps a possible way to obtaining 
an open fi tot < 1 universe from inflation. Guth [48-50] 
discusses the problem in detail, terming it the "young- 
ness paradox" . He speculates that this could explain why 
the SETI project has failed to discover extraterrestrial 
intelligence (since the preference for low cosmic age im- 
plies that we are overwhelmingly likely to be the first ob- 
servers in our Hubble volume), but finds it more plausible 
that it is merely a symptom of the probability calcula- 
tion technique not being the right one. A measure with 
i-ordering implies not only this, but also more extreme 
consequences such as that you personally are the first 
observer in your Hubble volume since thermalization, 



and that it is way more likely that you evolved on the 
first habitable planet to form, under circumstances where 
many apparent flukes accelerated your arrival. We con- 
clude that our living in a cool, leisurely universe means 
nothing other than that t-ordering is the wrong ordering. 

D. Why does not rule out eternal inflation but 
limits attractor behavior 

There is presumably some correct solution to the or- 
dering problem, and we have found that i-ordering is not 
the one. It is tempting to speculate that this is because it 
flaunts some important yet-to-be-established underlying 
physical principle, in which case the coolness problem is a 
helpful clue for identifying this principle. Here is a guess 
as to such a principle: Physical questions must be express- 
ible in terms of observables. This implies that it is a big 
no-no to condition on quantities that are not physically 
observable, like t or a. t-ordering violates this principle 
since, as described above, it corresponds to defining ref- 
erence objects such as, say, "protons at time t" . The 
global time t is a completely unobservable quantity, be- 
cause all we can measure in our thermalizcd region is 
time intervals, e.g., the time interval (t — itherm) since 
inflation ended here. Moreover, there are infinitely many 
alternative choices of time variable, and between causally 
disconnected regions of spacetime (such as two different 
thermalized pockets) , gauge ambiguities imply that there 
is no objective way of defining simultaneity. 

Let us now discuss the implications of the coolness 
problem for eternal inflation and attractor behavior. 
Since the same Fokker-Planck probability approach that 
elegantly predicts eternal inflation is ruled out by the 
coolness problem, is inflation really eternal in any mean- 
ingful sense? A standard definition of eternal inflation 
(e.g., [57]) is that "at any given time, part of space is in- 
flating, and the inflating volume increases over time ", yet 
this definition uses precisely the infamous time foliation 
by considering the situation "at any given time" . 

Inflation can nonetheless be eternal, in the following 
well-defined sense: A finite comoving volume can pro- 
duce an infinite physical volume and an infinite number 
of particles (and other reference objects). Equivalently, 
inflation can be eternal in the sense of producing infi- 
nite Level I multiverses according to the classification of 
[51]. This is readily seen in the following simple example, 
which is illustrated in Figure 4. Consider a single-field 
inflaton potential V(<p) with a peak (local maximum) at 
= 0, and a comoving spatial region where the inflaton 
potential varies extremely slowly in the x-direction: 

<f>(x, y, z) w ex, (93) 

where the constant e is so small that spatial gradients in 
<f> are negligible. If quantum diffusion is negligible, then 
inflation will eventually end everywhere in our comoving 
region except at the one point (f> — 0, but it clearly takes 
longer for <j> to roll off the peak and end inflation the 
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closer to <j> — it starts. In the slow-roll approximation 
and considering the very top of the peak where it can be 
approximated by an upside down parabola, one readily 
finds that the number of e-foldings required to roll down 
from <f> is given by (Appendix B.2) 

e^ocW" 1 , (94) 

where r\ = in 2 V " (0)/V(0) is the usual slow-roll parame- 
ter at the origin. The total thermalized volume produced 
by our comoving region is thus 

Voc J e 3N ^dVoc J e 3N ^dxoc J \x\ 3 ^ dx, 

(95) 

which will diverge near the origin and give an infinite vol- 
ume as long as r\ > — 3. This simply reexpresses the well- 
known result that slow-roll peaks are eternal inflation 
regions, but without involving the Fokker-Planck equa- 
tion or i-foliation. An alternative way of seeing this is to 
compute the probability distributions in equation (56). 
If the initial conditions give a uniform distribution for (f> 
around the peak, then since r\ < 0, equation (94) readily 
gives well-behaved probability distributions / lnf for AT to t , 
V cx e 3JVtot , a cx e^ ot and t ~ N to t/H n at thcrmalization, 
where H = m' 1 y / V(0)j3: 



P(N tot >N,) = e" w *, (96) 

P(V>V.) oc V? /3 , (97) 

P(a>a*) cx o2, (98) 

P(t>t,) w e r ' HoU . (99) 



Multiplying the corresponding differential probability 
distributions / mf by the volume factor w vo \ cx e 3Ntot in 
equation (56) gives the predicted probability distribution 
f p for the parameters of a random thermalized volume el- 
ement: 

P(N tot <N*) cx e ^ N ', (100) 

p(v<v.) cx v; 1+ " /3 , (ioi) 

P(a < a*) cx a 3+v , (102) 
P(t<U) °5 e (3+7?)if °*\ (103) 

Since 77 > — 1, these distributions are all unnormalizable, 
so the probability is zero that N tot , V, a or T is finite. 

The second type of eternal inflation region (the kind 
dominated by quantum diffusion) also survives even 
though the standard justification for it becomes invalid. 
The standard argument [32] is that quantum diffusion 
up the hill causes faster expansion that more than com- 
pensates for the fact that it would have been more likely 
to roll downward, but it is of course precisely such mea- 
sures favoring faster expansion that lead to the coolness 
problem. An alternative justification is to again consider 
the evolution of a single comoving spatial region, with- 
out applying any volume weighting until after thermal- 
ization. When quantum diffusion is included, cb(t) will be 



not the deterministic slow-roll trajectory, but a random 
walk. This is guaranteed to end after some finite time 
t on d when 4> leaves the slow-roll region of </>-space, after 
which the comoving region thermalizes with a volume 

Fcxe 3Wt °\ N tot = H{<f>{t))dt. (104) 

Jo 

Although both N and V are clearly finite for any trajec- 
tory (j){t), the volume expectation value (V) can easily be 
infinite. As a simple example, if the random variable N tot 
has an exponential distribution /jv(^Vtot) oc e ~ Ntot ^ N for 
some constant N (equation (96) is such an exponential 
distribution with N = —rj^ 1 > 1), then just as in the 
slow-roll peak case above, the probability that N to t is 
less than any finite value vanishes if TV > 1/3. Such an 
ensemble averaging is equivalent to what happens when 
averaging (integrating) over a finite comoving volume as 
in equation (95), since the quantum fluctuations are ef- 
fectively uncorrelated in causally disconnected patches, 
and infinitely many such patches develop. In other words, 
both of the "new inflation" mechanisms that can generate 
eternal inflation as per the ^-foliated Fokker-Planck for- 
malism can do so also in the sense of producing infinitely 
many reference objects in a finite comoving region. 

A different way of seeing this is to note that the ther- 
malization surfaces are perfectly well-defined physical hy- 
persurfaces in spacetime, so there are no ambiguities re- 
garding the question of whether they have finite or in- 
finite volume. Therefore the Fokker-Planck approach 
to computing this volume must give the same answer 
as the differential equation approach described above — 
indeed, equation (89) agrees with equation (103) with 
= (1 + t]/3)Hq. Ambiguities only arise if we ask 
whether one infinite volume is larger than another infi- 
nite volume. 

Although eternal inflation thus survives the cool- 
ness problem unscathed, the attractor dynamics giving 
testable prediction / p (p) independent of the initial con- 
ditions is severely challenged. There is of course nothing 
wrong with the Fokker-Planck equation per se: it can be 
consistently reexpressed in terms of any time variable one 
choses [8, 53, 62, 103-106], so the problem lies entirely 
with the limiting procedure used to extracting probabil- 
ities from it. However, modulo a tunneling caveat dis- 
cussed below in Section VUG, the Fokker-Planck equa- 
tion only exhibits this attractor dynamics when expressed 
in such a way that it rewards rapid expansion, since as 
described above, this was the mechanism by which the 
highest peak in the potential was able to dominate the 
asymptotic solution. To circumvent the coolness prob- 
lem, one must therefore avoid rewarding rapid expansion. 
There are two ways of doing this: 

1. Changing the reference objects to be comoving 
rather than physical volume elements, i.e., replac- 
ing the volume weight factor of equation (57) by 

Wvol(p) = 1- 
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2. Changing the time variable to something not re- 
warding faster expansion, like a or p. 

Either way, the attractor dynamics described above goes 
away, and the asymptotic distribution (if it exists) will 
depend on initial conditions (modulo tunneling issues to 
which we return in Section VII G). The former approach 
has been studied in, e.g., [8, 53, 62, 103-106], and can be 
argued to be unphysical since all reference objects that 
may be related to observers (particles, planets, etc.) are 
proportional to physical rather than comoving volume. 
Let us now consider the latter. Although changing the 
time variable to p or other observable quantities solves 
the coolness problem, we have seen that this leads to 
pocket-based orderings which we have already discussed, 
so let us turn instead to unobservable time variables such 
as a. 



E. a-ordering and others not favoring rapid 
expansion 

There are infinitely many possible choices of time vari- 
able, and different choices can give different predictions. 
[57, 58, 61]. As an interesting example giving qualita- 
tively different predictions from t-ordcring, let us con- 
sider ordering by the cosmic scale factor a. 

First of all, a-ordering is not afflicted by the cool- 
ness problem. This is readily seem from comparing 
equation (102) with equation (103) and noting that the 
exponential divergence is replaced by much milder one 
that is between quadratic and cubic. Since Tcmb = 
Tthermtttherm/d, repeating the derivation of Section VI F 
shows that equation (90) gets replaced by 

-PPcmb < T] = P[T thcim a thcim /a < T] 

= -P[athorm < aT/Tthorm] 

cx (aT/T therm ) 3 +" 

cx T z+r '. (105) 

The probability of observing Tcmb < 3K given 
than Tcmb < 4K therefore has the acceptable value 
(3/4) 3+,) w 0.5 rather than 10 _1 ° 56 as with t-ordering. 

As elaborated in, e.g., [57, 58, 61], replacing t by a 
as a time- variable reverses some qualitative conclusions. 
With the initial conditions of equation (93), the amount 
of thermalized volume produced up to a given time is 
readily found to be 



V thc 



,(S+r,)H t a 3+ v 



1 - 37T 1 



1 - 2,rj- 



(106) 



Since rj > — 1 at a the peak, the amount of thermalized 
volume produced thus grows steadily. Over time, the 
freshly thermalized volume comes from an increasingly 
tiny initial volume very close to the peak (reflected by 
the rj in the exponent), but this negative effect is more 
than offset by the huge inflationary expansion (giving the 
3 in the exponent). 



Equation (106) shows that when the inflaton poten- 
tial has multiple peaks, t-ordering and a-ordcring give 
radically different predictions. With a-ordering, taking 
a — > oo shows that the flattest peak (which has the 
largest [3 + 77]) will completely dominate the thermal- 
ized volume, and that peak height Ho is irrelevant. With 
t-ordering, on the other hand, taking t — > 00 shows that 
the peak with the largest value of (3 + r])H will com- 
pletely dominate the thermalized volume, i.e., that the 
most important factor is peak height H , modulated by 
a slight preference for flatness. By constructing new time 
variables that are functions of t, a and observables like 
H, one can readily construct a range of other predictions. 

In conclusion, global orderings based on a non- 
observable time variable are not all ruled out by obser- 
vation, but lack compelling physical motivation and, as 
a class, fail to make unambiguous predictions. Another 
argument against global ordering with "time" variables 
such as a, T, H and p is that they can change direction 
(when a spatial region starts to contract) and hence fail 
to smoothly cover the entire spacetime manifold. Most 
importantly, if we accept the observability principle of 
Section VII D that physical questions must be expressible 
in terms of observables, then this entire class of orderings 
is ruled out. 



F. Ordering by an observable and the pothole 
paradox 

Equation (62) shows that the number of reference ob- 
jects is proportional to the volume at thermalization, so if 
two initially equally large volumes roll off the same peak 
in V{<j>) in opposite directions, the one that requires more 
e-foldings to thermalize will produce more reference ob- 
jects in the end. As we will now see, this has an important 
implication for pocket-based orderings. 

Consider the number of reference objects that will ulti- 
mately form in a given comoving volume. Equation (62) 
shows that it is proportional to the physical volume at 
thermalization, i.e., that it grows during inflation and 
then stays constant. This implies that whatever time 
variable we use, most reference objects at a given "time" 
arc in regions that stopped inflating as recently as al- 
lowed by their existence. Above we saw that this ef- 
fect was rather mild with the time variable a but gave 
the catastrophic coolness problem with the time variable 
t. When conditioning not on a monotonically increas- 
ing time variable but on an observable quantity such as 
p, T or H (implying a pocket-based ordering), the ef- 
fect is that the longest possible inflation scenario will 
completely dominate the statistics, no matter how long 
it takes. This means that a flatter peak dominates over 
a shallower one. However, it also means that both will 
be outperformed if the inflaton potential has a pothole 
in the roll path where <f> will get stuck for an exponen- 
tially long time until it tunnels out [80] . The deeper the 
pothole, the better, so one could argue that almost all 
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observers can trace the origin of their pocket back to an 
exponentially rare tunneling event through a high peak 
in the potential. This event in turn would be overwhelm- 
ingly likely to have been preceded by a many other freak 
tunneling events. This is a disturbing scenario, since for 
any extremely long inflationary history there exists an 
even longer one, and and it is far from clear that it gives 
well-defined and observationally allowed cosmological pa- 
rameter predictions / p (p). Note that the exponentially 
small probability per unit time to tunnel out of a deep 
inflating minimum is irrelevant, since <fi is guaranteed to 
tunnel out eventually with probability unity. 

The "pothole paradox" above follows from strict ad- 
herence to weighting pockets by their number of ref- 
erence objects, even though each pocket contains in- 
finitely many 16 . In other words, it follows from assuming 
that some infinities are greater than others. This ar- 
gues against inter-pocket weightings 3, 4 and 5 from Sec- 
tion IV B 2, suggesting that that all countable infinities 
receive equal weight. 



G. Quantum tunneling and attractor behavior 

One area in need of further work is the effect of quan- 
tum tunneling in inflation. In some cases, for example 
ones involving the above-mentioned pothole paradox, this 
tunneling process may hold the key to understanding the 
qualitative predictions. 

As a specific example, consider the question of whether 
a-ordering gives attractor behavior or gives predictions 
that depend on initial conditions. As discussed above, a- 
ordering is equivalent to weighting by comoving volume 
V c , which is a conserved quantity. This means that, ig- 
noring quantum diffusion and tunneling, the distribution 
V c (<p, a) simply flows down into the valleys, never moving 
from one basin of attraction to another. In other words, 
initial conditions matter, since the probability of ending 
up in a given basin of attraction is simply the probability 
of starting out in that basin of attraction. Quantum tun- 
neling alters this conclusion, allowing transitions between 
basins either during inflation or after reheating, in the 
second inflation phase corresponding to a cosmological 
constant p\ > [114]. Modeling this as a Markov pro- 
cess with a matrix of jump probabilities per unit "time" 
a between the different basins, one can imagine conver- 
gence to an equilibrium state asa^oo where the basin 
weighting is independent of initial conditions. Further 
work should clarify whether such attractor behavior ac- 
tually occurs. For instance, although "downward" tun- 
neling to lower energy through the Coleman-de Lucia 
bubble nucleation mechanism is fairly well-understood, 



Note that global orderings lack this ambiguity, since there are 
only finitely many reference objects at any finite "time" (say t or 
a) as compared to an infinite number at finite "time" p, T, etc. 



the situation regarding "upward" tunneling is less clear. 
Also, there will be no tunneling out of a minimum with 
V{<p) = to first order, thereby spoiling attractor be- 
havior if there is more than one such minimum, but it is 
not clear that there is no tunneling out to second order. 



VIII. CONCLUSIONS 

If string theory or some other fundamental theory pre- 
dicts an effective inflaton potential V((f>), then it will 
clearly be worthwhile to compute the predicted cosmo- 
logical parameter probability distribution / p (p) to con- 
front it with observation. In this paper, we have at- 
tempted to strengthen the groundwork necessary for such 
a calculation in the context of classic slow-roll infla- 
tion, exploring how the predictions for the 8 parameters 
p = (f2tot; Pa, w , Q , n s , a, r, n t ) depend on both the infla- 
ton potential and the measure. We found that the results 
depend on a complex and interesting interplay between 
the two, except for the very simplest potentials with a 
single minimum that can only be rolled to in one way. 
It is a persistent myth that models where ensembles and 
anthropic constraints enter are untestable and hence un- 
scientific [112]. We have seen that, quite to the contrary, 
most such models that we have explored in the present 
paper are predictive enough to be eminently testable — 
indeed, so much so that they have already been ruled out 
by observational data! 

Since we found that the measure problem is currently 
the weakest link in the calculation of predictions, we will 
classify models below primarily by their potentials and 
only secondarily by their measure. This is intended not 
only to extend the shelf life of our discussion, but also to 
place the seemingly daunting measure problems in con- 
text by clarifying what they do and do not affect. 



A. How the inflaton potential affects the 
cosmological parameters 

The effective cosmological constant pa depends only on 
which basin of attraction <p rolls into, since it is deter- 
mined by the minimum in which <fi ends up. The five 
power spectral parameters (Q, n s , a, r, n t ) vary within 
the basin, since they depend on the direction from which 
<p rolls to the minimum — for one-dimensional potentials, 
there are only two possibilities (left and right), so they are 
determined uniquely by the half-basin (region between a 
potential maximum and minimum). The curvature pa- 
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rameter f2 t ot and total pocket volume oc e 3 tot depend 
on the starting point and hence crucially on the mea- 
sure. Finally, the dark energy equation of state w = — 1 
as elaborated below in Section VIII D. 

One of our key findings is a useful simplification: in all 
cases except m/j ~ mpi, the cosmological parameter pre- 
dictions are determined by extremum statistics alone. In 
other words, unless the characteristic scale A0 on which 
the inflaton potential V varies happens to be close to the 
Planck scale, the only aspect of V that matters obser- 
vationally is the statistical distribution of peaks and/or 
troughs. 

The rah 3> mp\ case predicts that 



1 r 

n s = 1--, 


(107) 
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a = ~32' 


(108) 
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nt = -g, 


(109) 


r w 0.15 ±0.02 


(110) 



and that the parameters (Q, p\) are determined by 
trough statistics, specifically by the joint distribution of 
V and its curvature at minima. Assuming that V = 
is not special in the fundamental theory, all that mat- 
ters is the curvature probability distribution at minima 
with V w 0, given for the one-dimensional case by the 
distribution for V"{<j>) where V'{4>) = and V{4>) w 0. 
The to/, <C Topi case predicts that 



n s ~ Oil, (111) 

a w 0, (112) 

n* w 0, (113) 

r w 0, (114) 



and that the parameters {Q,pA,n s ) are determined by 
peak statistics — for the one-dimensional case and V = 
not special, they depend mainly on the the joint proba- 
bility distribution for (V, 77) at the peaks. 

This means that aside from w and O to t, the number 
of precise and testable quantitative predictions is 3 when 
rah *C mpi and 4 when rah 3> TOpi, dropping to only 
one (n t = — r/8) when m/, ~ TOpi. In other words, in 
the 5-dimcnsional parameter space (Q,n s ,a,r,n t ), the 
predictions populate a 2-dimensional hypersurface when 
rah -C mpi, a 4-dimensional region when m/, <~ mpi and 
a 1-dimensional curve when mj, » mpi. In essence, the 
rah <C mpi limit simplifies since e — > and the rah ^> fh 
limit simplifies since {n s ,a,r) are all determined by N, 
which is in turn known to about 10%. 



B. How the measure affects the cosmological 
parameters 

Above we saw that except when nih ~ TOpi, the pre- 
dictions for the parameters (Q,n a ,a,r,nt,w,p\) could 



all be calculated from extremum statistics alone. In 
other cosmological examples where extremum statistics 
enter, such as the BBKS Gaussian-peaks formalism for 
galaxy clustering and the subsequent peak-patch formal- 
ism [113], the relevant extremum statistics can be ex- 
tracted unambiguously from the statistical properties of 
the function in question (in these examples, the cos- 
mic density field). In our case, however, these ex- 
tremum statistics are not necessarily those one would 
obtain by simply locating the appropriate extrema in 
V{<p) and tallying up the distribution. Rather, this 
"raw" distribution is modulated by the measure to give 
the above-mentioned extremum statistics that determine 
the parameter prediction. For example, for a homoge- 
neous Gaussian random field V (</>), the 4-dimcnsional 
vector [V{(f>),V ((/>), V" ((/>), V"'{(f>)\ has a simple multi- 
variate Gaussian distribution that can analytically be 
conditioned on V (((>') = and whatever else is appro- 
priate, but the measure (both initial conditions, ordering 
and the choice of reference objects) can multiply the re- 
sulting distribution of extremum properties by a factor 
that radically modifies it. Conversely, it is reassuring 
to note that modifying the extremum distribution is the 
only way in which the pesky measure problems we have 
discussed affect the cosmological parameter distribution 
unless rrih ~ mpi. 

Our simulated Measure A case essentially weighted 
all peaks/troughs equally (strictly speaking, troughs by 
the lengths of their basins of attraction and peaks by 
the lengths of the regions to which they flow), except 
for the rather straightforward conditioning on reference 
objects in Section VI. For other measures, the ex- 
tremum statistics can get highly weighted towards cer- 
tain peaks/troughs over others, either because of their 
intrinsic properties (like ^-ordering favoring height and 
a-ordering favoring flatness) or because of good neigh- 
bors (like a deep valley to tunnel from in the a-ordering 
case). We saw that t-ordcring is an extreme case where 
only the highest slow-roll peak matters, so that the prob- 
ability distribution is a 5-function corresponding to the 
highest peak (if m/, <C fh) or the trough that 4> can reach 
the fastest from the highest peak (if m,h S> fh). 

This author considers the inflationary measure prob- 
lem wide-open, so let us merely summarize our limited 
progress with it. We have argued that global orderings 
such as i-ordcring that reward rapid expansion are ob- 
servationally ruled out by the coolness problem (over- 
predicting the CMB temperature) . If we accept the no- 
tion that we may only condition on physically observable 
quantities {e.g., porT), then we are lead to pocket-based 
orderings that are specified by an intra-pocket ordering 
and an inter-pocket weighting. Radial intra-pocket or- 
dering provides a specific working recipe that may and 
may not be the correct one. Weighting all spatially in- 
finite pockets equally likewise provides a specific pre- 
scription that may and may not be right. It avoids the 
pothole paradox, but is computationally difficult to ap- 
ply because of the challenge of calculating the relative 
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frequencies with which different basins of attraction are 
rolled down into when quantum diffusion is taken into ac- 
count. This computation becomes still more complicated 
in recycling universe [114] and string landscape scenarios 
[16-29] where each thermalized region can trigger new 
inflation leading effectively to a Markov process of jumps 
between basins of attraction. 

Since there appears to be only a small-number of well- 
motivated measures and they make quite distinct param- 
eter predictions, it may be possible to settle the measure 
problem empirically by ruling out all but one of them 
observationally. 



D. Implications for dark energy 

We saw that in order for the standard inflationary cal- 
culation to be self-consistent (give the correct expansion 
rate H), the inflaton potential V(<p) must be defined to 
include any (^-independent vacuum energy contributions 
from other sectors of physics, including any bare cosmo- 
logical constant term that may be present in the Einstein 
field equations. The effective cosmological constant that 
we observe in our Hubble volume is therefore simply the 
height pa = V(<p a ) of the potential at the minimum that 
the inflaton has rolled down to here. Any messy inflaton 
potential V(<p) thus naturally predicts p\ ^ unless all 
its minima happen to have exactly the same height and 
this common height happens to be exactly zero. 



C. Implications for gravitational waves and 
spectral indices 

Great worldwide efforts are currently being devoted to 
measuring CMB polarization in the hope of measuring 
the primordial gravitational wave amplitude Q t via large- 
scale B-modes [115, 116]. A key theoretical question is 
therefore what amplitude Qt to expect, since this de- 
termines the prospects for upcoming experiments to dis- 
cover what they are looking for and influences their fund- 
ing and design. Our results are very encouraging in this 
regard. Models with rrih S> mpi predict r w 0.15, which 
(Figure 13) is close to current sensitivities and should be 
detectable already within five years by Planck + SDSS 
[85]. Models with rrih <C mpi predict unobservably small 
r-values, but appear to agree poorly with existing ob- 
servations: they fail to explain the near scale invariancc 
n s = 0.98 ± 0.02 that is observed [3, 86] by giving in- 
stead a wide distribution with — 1 < n s < 1 and they 
flounder on the smoothness problem by generally over- 
predicting Q. Finally, models with rrih ~ mpi typically 
predict r-values that are clearly within reach of a next 
generation CMB polarization satellite; Figure 13 shows 
typical values r ~ 0.03. 

The most imminent observational breakthrough may 
well be sharpening the current limits spectral index limits 
n s = 0.98 ± 0.02 to exclude the scale-invariant Harrison- 
Zel'dovich case n s = 1, and this will severely reduce the 
number of viable inflation models. We have seen that 
models with rrih S> mpi predict n s w 0.96 whereas mod- 
els with to/j <C mpi predict a broad distribution in the 
range — 1 < n s < 1 — i.e., they fail to explain why n s w 1 
and cannot give blue (n s > 1) spectra, the only exception 
being arguably unphysical measures that do not volume 
weight, like Measure B. 

Our rrih ?5 mpi predictions agree well with other ar- 
guments suggesting that classic single-field inflation pro- 
duces both gravitational waves and departures from scale 
invariance that will soon be observable [11, 117, 118]. 



1. Vanilla or not? 

If no new physics is added beyond this inflaton po- 
tential, a generic prediction is that this dark energy will 
exhibit "vanilla" properties indistinguishable from a cos- 
mological constant, i.e., w = —1, a time-independent 
density at late times (including the present) and no spa- 
tial fluctuations. The reason for this is that the effective 
inflaton mass scale m is huge by current standards 17 . 
This means that on a very rapid reheating timescale, <j) 
will settle down into its minimum and stay there, either 
promptly recollapsing space if the minimum is substan- 
tially below zero or giving p\(t) = V[(j)(t)] constant at 
late times. The huge mass scale m will not only ensure 
that the minimum is rapidly attained, but also that there 
are no observable dark energy fluctuations, since excita- 
tions in <f> are too massive to be excited by present-day 
energies. It is noteworthy that the substantial improve- 
ments (most recently [119]) in dark energy measurements 
since the first discovery [120, 121] have not established 
any departures from "vanilla" [122, 123]. However, im- 
proved dark energy observations are clearly crucial, since 
the fact that messy inflation may already explain pa ^ 
and its rough magnitude does not in any way preclude 
the existence of additional physics producing non- vanilla 
dark energy at late times, and there are interesting hints 
at modest statistical significance (e.g., [124]). As shown 
in Figure 3, improved dark energy observations also com- 
plement CMB polarization observations by measuring 
the exact same curve p(a) at a vastly lower density. 



If we Taylor expand a 1-dimensional potential around a V ~ 
minimum as V(4>) ~ m 2 <f> 2 /2, then equation (32) gives m = 
\/V" ~ mh 2 /m v . For the observationally favored case > 
mpi, equation (38) thus shows that m ~ Qmpi ~ 10~ 5 rapi ~ 
10 14 GeV. 
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2. Continuum or discretuum? 

Potentials V(<p) with infinitely many minima (such as 
Gaussian random fields as explored in Section V) can 
give a continuous probability distribution for p\ . In con- 
trast, the landscape picture emerging from string theory 
suggests a very large but yet finite number of minima, 
say n to t ~ 10 300 [16-29] 18 Suppose we find a fundamen- 
tal theory that has n sm minima which all reproduce our 
SU(3) x SU(2) x U(l) standard model of particle physics 
at low energies and differ only in their predictions for p\ . 
It has been argued that n am <C n to t, is., that only a 
tiny fraction of the minima reproduce our observed low- 
energy physics [125]. Our results show that the key for 
testable predictions is whether n sm 3> £ 4 Q 3 ~ 10 123 or 
not. If n sm 3> 10 123 , then there may be so many min- 
ima able to produce galaxies that pa will for all prac- 
tical purposes have a continuous probability distribution 
such as those plotted in this paper, i.e., the "discretuum" 
[16, 18, 19] of pA-values may be indistinguishable from 
a continuum. If n sm < 10 123 , on the other hand, then 
there may be only one or a handful of minima producing 
galaxies, raising the exciting possibility that a precision 
measurement of p\ would allow us to determine precisely 
which minimum we are in, predict the remaining decimal 
places of pa and bring closure to our quest for a funda- 
mental theory of physics. 



inflation models we explored typically overpredicted pa 
by a couple of orders of magnitude or more. This is di- 
rectly linked to the smoothness problem, where the same 
coupling between pa and Q multiplies the Q-distribution 
by Q 4 and overpredicts the CMB fluctuation level for 
most inflation models. In other words, the standard an- 
thropic explanation of the observed pA-valuc does not 
appear to work well unless we can solve the smoothness 
problem. 

E. Implications for spatial curvature 

Since almost all reference objects are in pockets of in- 
finite volume, all measures except arguably unphysical 
ones (that do not weight by reference object and hence 
do not favor infinite volumes over finite ones) predict that 
we should observe O tot w 1±10 -5 . A second implication 
is that the thermalized pocket that we inhabit is spatially 
infinite, i.e., that our Hubble volume is but one in an infi- 
nite "Level I multiverse" [51] with distant Dopplegangers 
and other Byzantine implications [51, 126]. The reader 
finding these implications disturbing can take heart in 
the fact that these inflation models remain quite falsifi- 
able, with Planck + SDSS forecast to shrink the currently 
allowed range fitot — 1.01 ± 0.02 [1, 2] by a factor of four 
to Afitot ~ 0.005 [85] and therefore potentially excluding 
this inflationary prediction. 



3. The smoothness problem 

It is striking that even though many people start 
foaming at the mouth when encountering anthropic pa- 
arguments, feeling that they render theories untestablc, 
we have found that requiring inflation models to predict 
a pA-value consistent with observation is in fact a very 
strong test that rules out broad classes of models. The 
key reason for this is that, as emphasized in [87, 98], re- 
quiring nonlinear structures to form strongly constrains 
the quantity pa/Q 3 whereas the selection effects acting 
on the fluctuation amplitude Q alone (with pa/Q 3 held 
fixed) appear to be substantially weaker. The suppres- 
sion of galaxy formation from increasing pa can therefore 
be largely offset by increasing Q, and we found that the 



F. Implications for quantum gravity 

If attractor dynamics makes the cosmological parame- 
ter predictions f p (p) independent of pre- inflationary ini- 
tial conditions, then inflation provides a cosmic censor- 
ship that erases all quantitative clues about preceding 
events at higher energy scales. We have argued that 
many measures providing such attractor behavior are ob- 
servationally ruled out by the coolness problem, which 
suggests that perhaps pre- inflationary conditions do mat- 
ter. If so, drawing conclusions about these conditions will 
not be easy, requiring both finding V{<p) and solving the 
ordering problem, but this nonetheless offers a glimmer 
of hope that precision measurements of cosmological pa- 
rameters may ultimately teach us something about quan- 
tum gravity. 



There may indeed be a hierarchy of minima in some sort of gen- 
eralized potential: 

1. Different ways in which extra dimensions can be compacti- 
fied 

2. Different discrete fluxes that stabilize the extra dimensions 
(this sublevel appears to be where the largest number of 
choices enter, perhaps 10 300 ). 

3. Once these two choices have been made, there may be a 
handful of different minima in the effective supergravity po- 
tential. 



G. Outlook 

In conclusion, we have studied the important but diffi- 
cult problem of how to compute testable cosmological pa- 
rameter predictions from complicated inflaton potentials. 
Although we have obtained a number of results summa- 
rized above, much work remains on multiple fronts - 
here are a few examples: 

1. Although many of our results apply also to the 
multidimensional case d > 1, others do not, and 
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a detailed study study of d > 1 would be of great 
interest — particularly d = 2 and very large dimen- 
sionalities like d > 10 2 . 



2. A better understanding of "upward" quantum tun- 
neling during inflation should settle the issue of 
if/when initial conditions matter. 

3. As concrete potentials V(4>) emerge from string the- 
ory or other fundamental approaches, it will be of 
great interest to repeat our calculations of / p (p) for 
these physically motivated potentials. 

4. Being the Achilles heel of the entire endeavor, the 
ordering problem is in grave need of further work. 
The urgency of this goes beyond cosmology, afflict- 
ing any ensemble theory with distinguishable infi- 
nite volumes. 19 

This will be challenging work. However, the avalanche of 
high-precision cosmological measurements is continuing 
unabated, so if the goal is to find an inflation model that 
predicts rather than postdicts the correct answer, then 
the time to rise to this challenge is now. 
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19 For instance, much recent work in the string theory landscape 
context has involved counting minima and quantifying their sta- 
tistical properties. However, it is far from clear with what sta- 
tistical weights these minima should be counted — equally (as 
done in most papers), by basin of attraction area or by some 
other measure related to, e.g., tunneling or inflation. Solving 
the measure problem may thus be required for string theory to 
qualify as a testable physical theory. 



APPENDIX A: COMPUTING THE 
OBSERVABLES FROM p(a) 

In this appendix, we review how the cosmological ob- 
servables are derived from the cosmic density history p(a) 
and derive the equations of Section II. This is mainly re- 
view material [8-15], presented here to clarify when and 
how most observables can be derived from p(a) alone 
without reference to V((p). 

Given the cosmic density history p(a), the expansion 
history a(t) is readily obtained by integrating the Fried- 
man equation 
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p : 



3m 2 



(Al) 



where the Hubble parameter H = d\na/dt = a/ a. Lu- 
minosity distances, angular diameter distances, and all 
other classical cosmological observables similarly follow 
from p(a), typically via integrals involving this function. 
We will now summarize the calculation of the observables 
Pa (a) , S]j(k), S] 1 (k) and the eight cosmological parame- 
ters from Table 1. 



1. Horizons and e- foldings 

Since the horizon radius is -ff -1 , the comoving horizon 
size a/H^ 1 = a oc l/a 2 p is constant on the dotted lines 
of slope —2 in Figure 3, increasing towards the upper 
right. Comoving wavelengths leave the horizon when- 
ever the comoving horizon size increases, i.e., whenever 
the curve p(a) is shallower than these dotted lines (when 
d\n p/d\na > —2, i.e., when d > so that the expansion 
is accelerating), and enter the horizon when the slope is 
steeper, d\np/ 'din a < —2. When two points lie on the 
same dotted diagonal, it therefore means that that the 
horizon volume at the two epochs is the same comoving 
spatial region. For instance, the comoving spatial vol- 
ume corresponding to our observable universe now at the 
present epoch (filled triangle in Figure 3) was the hori- 
zon volume during the inflationary epoch corresponding 
to the open triangle, and the perturbations that left the 
horizon then are those that we see entering our horizon 
now. 

Key epochs in Figure 3 are a sta rt (when inflation be- 
gan; a sta rt = if there was no beginning) a cx i t (open 
triangle, when our current Hubble volume left the hori- 
zon) , cik (when a fluctuation of wavenumber k leaves the 
horizon; open square or circle, say), a cn d (star, when in- 
flation ended), a rc hoat (cross, when reheating ended), a cq 
(filled square, matter-radiation equality) and ciq (filled 
triangle, now). The number of e- foldings before the end 
of inflation is 



(A2) 



Of particular importance is iVexit = N(a cx i t ), often re- 
ferred to as just N, which is simply In 10 times hori- 
zontal distance between the open circle and the star in 
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Figure 3. The figure illustrates that iV cx it depends not 
only on the behavior of p(a) during inflation, but also on 
when reheating ends. 20 The strong bounds on this quan- 
tity derived in [36, 37] can be readily read off geometri- 
cally from Figure 3. We must have p(a ex it) ^ 10 _12 mpi 
to avoid overpredicting gravitational waves and density 
fluctuations (as detailed below) and presumably need 
p(flend) > 1 GeV 4 for reheating to be able to produce pro- 
tons (c./. [127, 128]). This gives the extreme bounds are 
35 < iV 0X i t < 85, whereas most popular inflation models 
favor the narrower range -/V cxit ~ 55 ± 5 [36, 37]. 



2. Spatial curvature 

The quantity |fitot — 1| -1 ^ 2 is simply the ratio of the 
curvature radius (oc a) to the horizon radius The 
spatial curvature parameter Qtot therefore evolves as 
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i.e., f2tot stays constant on the dotted lines of slope —2 
in Figure 3 and increases towards the upper right. This 
means that fitot is identical on horizon exit and subse- 
quent horizon entry, so that our currently observed f2 tot 
equals fitot(aexit)- We will make the standard assumption 
that, if inflation had a beginning, things were a mess at 
that time with fluctuations of order unity on the horizon 
scale, i.e., f2tot(astart) ~ 1< The inflationary prediction 
is therefore 



where iV bo fore = ln(a st art/aexit) is the number of c- 
foldings of inflation before our comoving Hubble volume 
exited the horizon (open triangle in Figure 3) . In typical 
slow-roll inflation models like the one plotted in Figure 3, 
the p(a)-curve is nearly horizontal and so equation (A3) 
is totally dominated by the second term, giving simply 
l^tot — 1| ~ e _2Arboforo . To match the current observa- 
tional constraint |O tot — 1] < 0.02 [1, 2] therefore requires 
-^before ^ 2.3. If the true curvature were exactly zero, 
density fluctuations on the horizon scale would propagate 
into our measurement and give |f2tot — 1| ~ Q ~ 10~ 5 , 
corresponding to -/Vbefore ~ 5.8, so the prediction for the 
observed curvature is 
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FIG. 17: Zoom of the lower right corner of Figure 3: 1<t constraints 
on the densities of matter and dark energy computed from SN la 
[119], CMB and galaxy clustering data as described in [122]. 



and Abefore-values in the range 6 < -^before < 00 & re 
observationally indistinguishable. 



3. Dark energy 

The dark energy density is, in practice, defined as the 
part of the observed density that we do not understand. 
The density contributions that we do understand at least 
at some level are those from photons, matter and curva- 
ture, py, p m and pk: 



p 7 (a) 
Pm(a) 
Pk(a) 
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where T = -^T is the CMB temperature and C(3) 
1.202. The dark energy density is thus the residual 



As discussed in [36], if <p oscillates around a generic (parabola- 
shaped) minimum in the potential during reheating, the loga- 
rithmic slope dlnp/dlna = —3 just like during matter domina- 
tion, so ending reheating at higher density increases N Ex j t . If 
the second derivative of the potential vanishes in the oscillation 
direction(s), however, as for V{<j>) oc </> 4 , then dlnp/dlna = —4 
like during radiation domination and N ex - lt becomes independent 
of when reheating ends. 



PA (a) = p(a) - p 7 mk(a), (A9) 

where p 7m k = p 7 + p m + Pk- Since all parameters enter- 
ing into p 7 mk(a) are now fairly well measured (f2tot, £b, 
£ c , as in Table 1, H , T ) [1-3], equation (A9) deter- 
mines pa (a) directly from the behavior of p(a) around 
the present epoch. If <fi settles into a stable potential 



38 



minimum at (p a long before the present epoch, then we 
observe p\ = V(<p Q ) independently of time, som = — 1. 

Dark energy is often discussed in the literature as 
something separate from inflation, so our merging of the 
two in this paper requires an explanation. 

First of all, note that if General Relativity has a true 
cosmological constant /3^ R hard-wired in, quite separate 
from the inflaton-only potential V((p), then we will only 
obtain the correct predictions from inflation if we absorb 
this cosmological constant into the potential and insert 
V(cp) = V(cp) + pa into the equations below. Inserting 
V{4>) would give the wrong answer, since the Hubble 
parameter plays a key roll during inflation and depends 
on the total density via equation (Al), including p^. 

Second, note that if there is some other mechanism 
generating a "dark energy" density contribution without 
having any connection to the physics of inflation, it must 
be absorbed into V(tp) for the same reason. Moreover, 
for scalar-field models of dark energy, at least ones with 
standard kinetic terms, it is mathematically straightfor- 
ward to perform this inclusion, cither by modifying the 
dynamics of the inflaton <fi at very low energies or by in- 
cluding the dark energy fields as additional components 
in the vector (p. 

Third, generic potentials can have many minima, typi- 
cally with V((p) 7^ 0, so they can readily produce nonzero 
dark energy regardless of whether there is also a non- 
inflationary source of dark energy. If inflation is related 
to the string landscape, say, then inflation-based dark en- 
ergy may not be optional: if there are say 10 300 minima 
[18-29, 31] and > £~ 4 <2~ 3 ~ 10 123 of them reproduce 
our low-energy physics, then it is likely that many min- 
ima have pA-values in the range |pa| ^ £, 4 Q 3 permitting 
galaxy formation. 



larger k give larger at (open square and circle). During 
the observable part (N <~ 55) of slow-roll inflation, H 
rarely drops by more than about a percent per e-folding, 
so for cosmologically observable scales the approxima- 
tion H(ak) ~ H(a cx it) is quite accurate, giving simply 
as, = aok/H(a ex n) and thus 

This has a simple geometric interpretation in terms of 
Figure 3. Since 5^ oc p 1 / 2 and k oc ak, the curve Sj^(k) 
in a log-log plot is simply the p(a) curve from Figure 3 
compressed by a factor of two vertically. In other words, 
Stfik) and p(a) are the exact same curve, merely with 
different axis labeling. 

This means that tensor spectral index 

^ _ din Si 2 _ riln<% 2 = dlnp 
dink rflna dlna 

is simply the slope of the p(a)-curve in Figure 3, evalu- 
ated at a cx ; t (open triangle). 

5. Density fluctuations 

None of the three above-mentioned observables f^, 
PA (a) and 5]j(k) depended explicitly on the physics of 
inflation, being computable directly from the curve p(a). 
We will now see that this is in many cases true even for 
the fourth observable, the density fluctuations fe(fc). 

a. Computing p(a) and &(&) from V(4>) 



4. Gravitational waves 

During inflation, quantum fluctuations in any mass- 
less field ip will produce effectively classical fluctuations 
with power spectrum P^(k) = (H/2ir) 2 at the epoch 
when the scale k exits the horizon [11]. The modes of the 
spacetime metric corresponding to gravitational waves 
are examples of such massless fields, and therefore pro- 
duce current gravitational wave fluctuations [11] 
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Note that this has nothing to do with fluctuations in the 
inflaton field, depending on the cosmic density history 
p(a) alone. 

The exit scale a^ is computed by solving the equa- 
tion akH(ak) = aofc for a*,, where H = (p/3) 1 / 2 /m as 
per equation (Al). In other words, is determined 
by following the dotted diagonal in Figure 3 from the 
point (a, p) = (aok/Ho, p(ao)) up to the left until it inter- 
sects the inflation curve. The current horizon wavenum- 
bcr k = Hq thus gives ak — a xit (open triangle) and 



The cosmic density is given by [11] 

1, .. 



p = V{<p) + -\<pY 



Pjmk 



(A13) 



and the evolution of the inflation vector is governed by 
[11] 



+ 3i?0 + W(0) = 0. 



(A14) 



(For a review of inflation generalizations with non- 
standard kinetic terms, see [11].) We ignore the p 7m k- 
term in equation (A13) during inflation since it rapidly 
becomes negligible, having a steeper slope (< —2) than 
that of p (> -2) in Figure 3. 

The slow-roll approximation (SRA) is that the </>-term 
in equation (A13) and the 0-term in equation (A14) are 
negligible, which gives 



p = V, 



<P=- 
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(A16) 
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Using equation (Al) and the identity d/dt = 
Hd/dlna = —Hd/dN, this becomes 



In other words, the slow- rolling inflaton vector evolves 
just like a particle moving in an extremely viscous 
medium, so that the velocity is proportional to the force, 
and it always rolls in the direction of the gradient, i.e., 
takes the steepest descent. The only differences are that 
the time variable is — N rather than t and that the po- 
tential is In V rather than V. For the 1-dimensional case, 
equation (A17) is usually integrated as 



N = In- 
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If valid, the SRA remains valid as long as the slow-roll 
parameters e«l and \rj\ <§; l 21 , where 



e=^|VlnV| s 
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For the multidimensional case, V" denotes the second 
derivative in the direction in which the field is rolling, i.e., 
in the direction of the gradient V<f>. For the 1-dimensional 
case, the third slow-roll parameter is defined as 
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For the (i-dimensional case (where the vector cp has d 
components), the power spectrum of density fluctuations 
S H (k) is given by [10, 11, 129-131] 
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where the right hand side is evaluated at the exit epoch 
afe as usual, and N(<pi) denotes the number of c-foldings 
for the field to roll from <f> down to some fixed reference 
density according to equation (A17). Heuristically, equa- 
tion (A21) is easy to understand: quantum fluctuations 
in cp are of order Sep ~ H <~ p 1 / 2 /fa, and if such a fluc- 
tuation in some particular spatial region delays the roll 
down to our current epoch by an amount 5N ~ ^^S(f>, 
this will appear an an overdensity of order SN. Fluctu- 
ations from the different components <j>i of the inflaton 



21 For the multidimensional case, \rf\ C 1 is strictly speaking not 
a sufficient condition. The second derivative <p can get large 
either because the inflaton starts rolling faster (large \n\ so that 
| V In V\ grows) or because it changes direction (a large eigenvalue 
of the Hessian matrix (InV), ij corresponding to an eigenvector 
misaligned with the gradient) . Generically, such eigenvectors will 
have a component in the gradient direction as well and cause a 
large |?7|, so it requires some fine-tuning to fool the \n\ and e tests. 



vector cp add in quadrature to give the total variance 
Sfiik) 2 ■ Choosing our coordinates in </>-space so that the 
field is momentarily rolling in the 0i-direction, we have 
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In most studied cases, the walls of the multidimensional 
gorge in which the inflaton slowly rolls are much steeper 
than the roll direction, rendering the inequality of equa- 
tion (A22) close to an equality [10, 11, 129-131]. 



b. Computing 8h (k) from p(a) 

Let us now eliminate V{cp) and rewrite the above quan- 
tities in terms of the p(a) alone for the one-dimensional 
case. Letting primes ' denote d/dN = —9/9 In a, we 
obtain 

(A23) 
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again evaluated at the exit epoch cifc. This gives 
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In deriving these and related relations, the following 
identities arc useful (for brevity, we here — and only 
here — we use the shorthand notation ' = ^ and and 
r = In p = In V ) : 



d N 
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-m 2 rd^ 



8% = rh 4 (r 2 d 2 +rrd^), 
8 3 N = ~fh 6 [f 3 9| + 3f 2 fd 2 + (r 



- -9 . 9 ■ • ■ 

rr + r r 



(A30) 
(A31) 
(A32) 
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etc. The derivatives of a curve (In a, In p) = (—N,r) 
therefore satisfy 



r 

r" 
r"' 
Jin 



-2-2 

—to r , 
2m 4 f 2 f , 



-m 6 f 2 (2f 2 + ff) 

-2m 8 f 2 (4f 3 + 7rr"f + f 2 'r') 



r = 



to 



1 r" 
1 



m 3 (-r' 3/2 ) 



2r' 



(A33) 
(A34) 
(A35) 
(A36) 

(A37) 
(A38) 
(A39) 



and the cosmological predictions of equations (37)- (43). 
2. Inflation near a maximum 



Taylor expanding V((f>) to second order around a max- 
imum O gives 



U(0) = V(<h) 



1 + 



r/o ( - O 



(B8) 



where rjo is the second slow-roll parameter of equa- 
tion (A19) evaluated at the maximum. Generically, 
r/o < 0. Using equation (A19), this gives 



APPENDIX B: INFLATION NEAR PEAKS AND 
TROUGHS 

Since the cosmological parameter predictions for in- 
flation with rrih 3> to or rrih <C to depend only on the 
behavior of the inflaton potential near its minima or max- 
ima, respectively, let us review the analytic solutions for 
these two important cases and compute the predicted 
cosmological parameter probability distributions. 



e = 



2m 2 



Vo ( - <h 



v = vo 



1 , m ( </>-<A) 

2 V m 



2 \ m 

-2 



(B9) 
(BIO) 



In the early stages of <j> rolling off the peak, while the term 
in square brackets remains near unity, equation (A18) 
gives 



Inflation near a minimum 



(<j> - 4> ) oc e 



VoN 



(Bll) 



Taylor expanding V(<p) to second order around a min- 
imum <j>o gives 



Vtf) 



V(<f> ) + -V"(cb )(<j>- 



,) 2 



(Bl) 



where, generically, V"{<j>) > 0. Since we saw in Section VI 
that only those minima where \p\\ = \V((f>o)\ was neg- 
ligibly small contributed significantly to the parameter 
probability distribution, we need only consider the sim- 
ple case V((f>o) — 0, i.e., the standard parabolic potential 
case [32]. Using equations (A19), this gives 



e = J] = 



2m 2 



(0-<M 2 



so inflation ends with e = 77 = 1 when 
to 



fond 



k) ± V2-. 



m. 



(B2) 

has rolled down 
(B3) 



Equation (A18) now gives equation (37), which in turn 
(via equations 11-20) gives the well-known results 



N = 



M 2 



477J 2 



(f> = 0o ± fh y / 2(2N + 1) 
1 

e = 77 = 

6 = 0, 



2N + r 



(B4) 
(B5) 
(B6) 
(B7) 



so at the observable epoch TV ~ 55, the field is expo- 
nentially close to the peak unless rjo happens to lie in 
the small interval —0.02 < 770 < (slow-roll requires 
— 1 < rjo < 0). According to the 2nd order expansion of 
equation (B8), inflation ends only when \<f> — cj> \ > to. 
However, the scaling relations given in Section V show 
that inflation ends much sooner because of higher or- 
der terms in the Taylor expansion, usually because the 
cubic term causes inflation to end when \rj\ = 1 with 
1 4>cnd — ^0 1 ~ to/j 3 /to 2 , which is borne out by the rrih <C to 
simulations in Section V. For the observable regime 
N <~ 55, we thus have 



l^-^ol 

V(4>) 

U(0) 



V 

6 



mh e noN 



Va m h 6 c 2 noN _ q 

TO 6 



= TO 



V 2 



rioN 



0, 



(B12) 
(B13) 

(B14) 
(B15) 
(B16) 



and the cosmological predictions of equations (48)- (54). 
In equation (B16), we used the scaling relation V'" /V ~ 
to/i 3 which follows from equation (32). 
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